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FOREWORD 


The results presented in this report were the accomplishments of a 
research and development project performed for the Vibration and Acoustics 
Branch of the Goerge C. Marshall Space Flight Center, National Aeronautics 
and Space Administration under Contract No. NAS8-21403. This work has been 
accomplished by the Vibration and Acoustics Group, Structural Engineering 
Branch, Chrysler Corporation Space Division, Huntsville Operations from 1 
July 1968 to 31 July 1969. Eleven Monthly Progress Reports have been sub- 
mitted for the past months presenting derivations of equations and analyses 
of results. The final results of this project are presented in Two Volumes: 
This is Volume One, Theoretical Analyses, and Volume Two is a Users' Manual 
for Computer Program RANDOM. 
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ABSTRACT 


Formulations were derived and computer programs were written to 
calculate the random vibrational responses of rectangular cylindrical 
panels cross-reinforced with ribs and stringers subjected to the fluctuating 
pressure enviornments . Three boundary conditions are considered: four 

edges simply supported; four edges clamped; two opposite edges simply sup- 
ported and the other two clamped. Special cases of complete cylinders and 
flat panels are also included. The computer programs can be selected to run 
on any one or all three of the boundary conditions. Either the spectral 
density or the one-third -octave level of the excitation pressure may be 
input in any discrete frequencies. The computer programs will apply when 
either the complete panel or a portion of it is subjected to the excitation 
pressure, The frequency equations are incorporated in the computer programs. 
Input data include only the geometric dimensions and material properties of 
the panel. Formulations are according to the normal mode approach. The 
responses calculated are the acceleration, the displacement, and the stress 
spectral densities. Mean-square and root-mean-square values are calculated by 
numerical integration. All spectral densities are tabulated and plotted. 

New expressions for the joint acceptance of all mode combinations for dif- 
ferent correlation functions are derived. Both the main and the cross terms 
are taken into account to obtain the responses. Up to 625 terras are summed 
to obtain the spectral density at one data point. More than 1000 data points 
are calculated for each spectrum. Both the local responses at any point and 
the average responses over the complete panel are calculated. The program 
can be used to investigate the contribution of main terms and cross terms. 
Comparison of calculated results with test data shows good agreement. 
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I . INTRODUCTION 


The purpose of this project was to develop computer programs to calculate 
the random vibrational responses of rectangular cylindrical shell panels cross- 
reinforced with ribs and stringers. The boundary conditions considered are four 
edges simply-supported, four edges clamped, and two opposite edges simply-sup- 
ported while other two clamped. The special cases of complete cylinders and 
flat panels are also included. 

A total of more than ten programs have been developed. The main program 
RANDOM contained 'the formulations for the three boundary conditions and will 
calculate the responses with a single loading of the input data. It can also 
be selected to calculate the responses for any one of the three boundary 
conditions. For accurate results, and large frequency range, programs were 
written for each boundary condition to perform special investigations. 

The one-third octave spectrum of the excitation pressure is read in as 
input in a sequence of any discrete frequencies. Hence, any excitation pres- 
sure spectrum of any shape can be used as input. This improves the simulation 
of the excitation pressure field. The excitation spectrum is converted into 
pounds-per-square-inch squared per Hertz. Excitation pressure for each data 
point frequency is obtained by interpolation. The excitation pressure spec- 
trum is plotted out both in decibel scale and in (psi)2 per Hertz. 

The formulations are according to the normal mode approach. Both Alan 
Powell's El] joint acceptance and Y. K. Lin's [2] cross spectral density of the 
generalized force are used in the formulations. The relation of these two 
quantities is derived. The advantages and disadvantages of each are discussed. 
The analytical expressions of these two quantities for all mode combinations 
for two correlation functions are derived. These are believed to be new. One 
correlation function is exponentially decaying with separation distance and 
frequency while the other is a cosine function with exponentially decaying 
amplitude. Separate computer programs are written to study the joint acceptance 
and the cross spectral density of the generalized force. It is found that the 
curves of the joint acceptance and the cross spectral density of the generalized 
force are very regular for the exponentially decaying correlation function. The 
behavior of these quantities for the sinusoidal and exponentially decaying 
correlation function is very irregular. 

Contributions of both the main terms and the cross terms are summed to 
obtain the responses. The one n-th octave bandwidth is used for frequency incre- 
ment to save computer time and yield smooth spectral density plots of the 
responses. The frequency range for the spectrum is 5000 Hertz or more. Up to 
625 terms are summed to give the response spectral density at each data point. 
More than 1000 data points can be calculated for each response spectral density 
plot. The responses are calculated as the displacement, the acceleration, and 
the stress spectral densities. Mean-square and root— mean-square values are 
calculated by numerically integrating the area under the spectral density curve. 
The response spectral densities are tabulated and plotted in graphical forms with 
the root-mean-square value printed at the top of the plot. The programs will 
apply when either the complete panel or a portion is exposed to the excitation 



pressure. Both the local responses at any point and the average responses over 
the whole panel can be calculated by the programs . 

In one of the programs, the acceleration spectral density is expressed 
in decibels referenced gravity acceleration and the vibro-acoustic transfer 
function is calculated as the acceleration spectral density minus the excitation. 
This is very useful for the investigation of the transfer function of aerospace 
vehicle structures. 

Separate programs are written to investigate the contribution of the 
cross terms to the total response. It is found that although the cross terms 
do not contribute very much to the mean-square response, they do affect the 
shape of the response spectral density plot to a certain degree. 

The natural frequencies of the panel are calculated in the programs. 

These frequency equations are newly derived or modification of those available 
in the literature. For panels of uniform thickness, these equations are the 
same as those in the literature. For cylindrical shell panels cross-reinforced 
with stiffeners, no frequency equations can be found in the literature that can 
be advantageously used in the programs. These newly derived equations though 
approximate in nature, yet yield reasonable results. The frequency equations 
account for the boundary conditions, the rigidity of the stiffeners, and the 
curvature of the shell. They are not very complicated so they can be incor- 
porated into the computer programs without requiring a large amount of computer 
time. 


A program is written to calculate the total number of natural frequencies 
up to a certain range and the modal density. This program is useful in the 
investigation of the dynamic characteristics of the structure. 

By utilization of the developed computer programs, investigations on the 
effect of boundary conditions on the responses, on the contribution of each 
mode and the cross terms to the responses are performed. It is found that the 
more rigid the boundary condition, the less the displacement spectral density 
and the larger the acceleration spectral density at resonance. The root-mean- 
square displacements for the three boundary conditions are not very much different. 
Estimation of the responses are made by the spectral density at fundamental mode. 
It is found that the fundamental mode contributes up to 50% of the mean— square 
responses, and the second mode contributes more to the acceleration than the 
displacement responses. 

By assuming that the continuous structure vibrates at its fundamental 
mode, simple formulas are derived to calculate the response spectral density 
and the mean-square responses. It is found that good approximation of the spec- 
tral density at fundamental mode and the mean-square responses can be obtained by 
intelligent use of these formulas. 

Comparison of the computed results with experimental data from projects, the 
Chrysler Huntsville Operations is presently conducting for Marshall Space Flight 
Center, shows good agreement. 
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Since all the derivations of the equations used were reported in the 
Monthly Progress Reports, only the resulting equations are presented in 
this Final Report. For detail derivation, the reader will be referred to 
the applicable Monthly Progress Reports. 
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II . FORMULATIONS 


2. 1 Frequency Equations 

2.1.1 Four Edges Simply-Supported Rectangular Cylindrical Shell 
Panels Cross-Reinforced With Stiffeners 


The undamped natural frequencies of four-edges simply supported 
rectangular shell panels cross-reinforced with ribs and stringers may be ex- 
pressed as 


01 


jk 



+ 



( 2 . 1 - 1 ) 


2.1.2 Four Edges Clamped Rectangular Cylindrical Shell Panels Cross- 
Reinforced with Stiffeners 

The frequency equations for cylindrical shell panels with 
edge conditions other than simply supported are involved and complicated. 
However, reasonable accurate expressions may be given as follows: 


ai 
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2.1.3 Two Opposite Edges Simply-Supported and the Other Two Clamped 
Rectangular Cylindrical Shell Panels Cross-Reinforced With 
Stiffeners 

\ 

The frequency equations for this case are as follows : 



5 




j — 2,3, • • . 


where 

“jk ~ undamped natural frequency in radians per second 


(2.1-3) 
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The smeared-out mass per unit area is 

M = ph+p'h' (2.1-5) 

E = Young's modulus of panel skin 

E’ = Young's modulus of stiffeners 

h = Thickness of panel skin 

h' = Smeared-out thickness of stiffeners 

v = Poisson's ration of panel skin 

a = Radius of shell 

= Spacing of width-direction stiffeners 
b-|_ = Spacing of length-direction stiffeners 

b = Circumferential width of panel 

l = Axial length of panel 

1 2 = Moment of inertia of one length-direction stiffener with respect 

to neutral axis 

I2 = Moment of inertia of one width direction stiffener with respect 
to neutral axis 

p = 'Mass density of panel skin 

p ' = Mass density of stiffeners 
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2.1.4 Flat Panels 


The flat panel is the limit of a cylindrical shell panel when 
the radius approaches infinite. All the frequency equations for the shell 
panels reduce to the corresponding equations of the flat panels, when the 
radius "a" approaches infinite. 

2.1.5 Complete Thin-Walled Cylinders 

It can be shown that the frequency equations of the two- 
side simply-supported and the other two clamped rectangular cylindrical shell 
panels reduces to those of complete thin-walled cylinders with both ends 
simply-supported when the circumferential width "b" is equal to ra. Similarly 
the frequency equations for complete thin-walled cylinders with both ends 
fixed are identical with those of four-side clamped rectangular cylindrical 
shell panels with circumferential width equal to 7ra. Hence, the upper limit 
of the circumferential width for the validity of the frequency equations is 

b _< ira 

When b equal to ira, the difference between the complete cylinders and the 
shell panels lies in that the index k in the frequency equations gives the 
number of waves around the circumference of the complete cylinder but the 
number of half-waves in the shell panel. The index j gives the number of 
half-waves in both the complete cylinder and the shell panel. 

2.1.6 Cylinders and Shell Panels of Uniform Thickness 

When the moment of inertia of the stiffeners, [Ip and Ig in 
Equation (2.1-4) ] approach zero or when the spacings of the stiffeners [ap and 
bp] in Equation (2.1-4) approach infinite, the cylinders and the shell panels 
will be of uniform thickness. Equation (2.1-4) shows this is actually the case. 
Therefore all the frequency equations apply to both structures of uniform thickness 
as well as reinforced by ribs and stringers . 

2.2 The Normal Modes 

The normal mode shapes of rectangular cylindrical shell panels cross- 
reinforced with ribs and stringers subjected to different boundary conditions 
may be reasonably represented as 

F jk (r) = X j (x)Y k (y) (2.2-1) 


The functions 
sections . 


X j 


and Y, 


for each boundary condition are given in the following 
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2.2.1 Four Edges Simply-Supported Rectangular Cylindrical Shell 
Panels Cross-Reinforced With Stiffeners 

For the four edge simply-supported rectangular cylindrical shell 
panel cross— reinforced with ribs and stringers, the Xj and Y^ functions may 
be reasonably represented by the normal modes of four edges simply supported 
rectangular flat panel of uniform thickness, namely: 


= sin 

j irx 
£ 

(2. 2- 2a) 

= sin 

k7f y 

b 

(2. 2- 2b) 


2.2.2 Four Edges Clamped Rectangular Cylindrical Shell Panels Cross- 
Reinforced With Stiffeners 

. For the four edge clamped rectangular cylindrical shell panel cross- 
remforced with ribs and stringers , the Xj and functions for the normal modes 
may be reasonably represented by 


1! 

X-^(x) = cosh 

1.5056irx 

a 

_ 1.5056ifx 

COS £ 

(2 . 2-3a) 

- 

0.9825 ^sinh 

1.5056t rx 
l 

sin 1.5056irx\ 
l / , 


Yi = 

Y^(y) = cosh 

1.50567ry 

b 

cos liLMiy 
b 

(2.2-3b) 

- 

0.9825 ^sinh 

1.5 0 5 67ry 
b 

1.5056iry\ 
sxn ^ 


li 

Xj (x) = cosh 

(j + 1/2)t7x (i + l/2)rx 

£ ~ COS £ 

(2.2-3c) 

_ 

sinh (l + 1/2)ttx + 
£ 

3 = 2, 3, ... 

(.1 + l/2)irx 

i 


Y k “ 

Y k(y) = cosh • 

(k + l/2)ny 

b 

- cos ms. 

(2.2-3d) 

- 

sinh + gin 

b 

(k + 1/2) tty 
b 



k = 2, 3, . 
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2.2.3 Two Opposite Edges Simply-Supported While Other Two Clamped 
Rectangular Cylindrical Shell Panels Cross-Reinforced With 
Stiffeners 

For two opposite edges (perpendicular to x-axis) simply- 
supported and other two clamped rectangular cylindrical shell panels cross- 
reinforced with ribs and stringers, the Xj function is given by Equation 
(2.2-2a) while the functions is given by Equations (2.2-3b) and (2.2-3d). 

2.3 Spatial Correlation of the Excitation Pressure Field 

The spatial correlation coefficient of a time stationary excitation 
pressure field at any two points is defined as 


where 


-T- -> < p(r 1 )p(r 2 )> 

p(r;L ,r 2 ) = 2 


(2.3-1) 




P ( r ^) jP (r 2 ) 


< P(r 1 ),p(r 2 )> 


2 

?r 


pressure at the two points r-^ and r 2 
temporal average of the pressure product 


/ 


V <p 2 ( r 1 )> <p 2 (r 2 )> 


reference mean-square pressure 


For space stationary excitation pressure fields, the spatial correlation 
coefficient will be a function of the separation distance 


where 



<P (r 1 )p(r 2 )> 




(2.3-2) 


For the stationary both in time and space pressure field, the mean-square 
pressure of the field will be used as the reference pressure. 

For the fluctuating pressure field encountered in aerospace vehicle, both 
theory and experiments indicate the correlation coefficient may be reasonably 
represented by the following expressions 
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= expC-A-jKC) 

(2.3-3) 

(5) 

= exp(-A 2 K£) cos (Kg) 

(2.3-4) 

5 

= separation distance 


A 1’ A 2 

= decay constants 


K 

- — = wave number 

c 


0) 

= circular frequency 


c 

= speed of the pressure wave 



Equation (2.3-3) expresses the correlation coefficient decreasing exponentially 
as the separation distance and wave number increase. This is appropriate to 
be used when the pressure field is a boundary layer field. Equation (2) ex- 
presses the correlation coefficient as a cosine function of the product of the 
wave number and separation distance with amplitude decreasing exponentially. 
This will apply to progressive pressure fields. 

Typical curves for equation (2.3-3) are shorn in Figures 43 and 44. In 
Figure 43: 


5 


hi 


C 


In Figure 44: 


A 


1 




100 Hertz = 200ir radian/second 
13500 inches /second 


0.5 


100 Hertz = 200?r rad/sec 


c = 13500 in/sec 
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Typical curves for Equation (2.3-4) are shown, in Figures 45 and 46. In 
Figure 45 : 


A 2 = 5 

w = 100 Hertz - 200t rad/sec 

c = 13500 in. /sec. 

In Figure 46 : 

A 2 = 0.5 

«J = 100 Hertz = 200 tt rad/sec 

c = 13500 in. /sec. 


2.4 The Joint Acceptance and the Cross Spectral Density of the Generalized 
Force 


New expressions for the joint acceptance squared and the cross spectral 
density of the generalized force of rectangular panels subjected to the 
excitation of fluctuating pressure environments are derived in the following. 
In this derivation, the analytical expression of the cross spectral density 
function of pressure field conforms quite accurately to the actual environ- 
ments of the aerospace vehicles. Analytical expressions for the joint accept- 
ance squared and the cross spectral density of the generalized force as a 
function of frequency are derived for all possible mode combinations. This is 
believed to be new. 

2.4.1 Cross Spectral Density of Generalized Force, Sinusoidal Decaying 
Correlation 

The cross spectral density of the generalized force of the pres- 
sure field is defined as [see Reference 2, Chapter 7, Section 7-2, Page 211, 
Equation (7-39)]: 


jkmn 

= / / V ( ^l’ ^2’ w) F jlA> Wr 2 ) & 2 

(2.4-1) 

F 

mn 

= normal modes of the structure 



For space homogeneous pressure field, the cross spectral density is equal to the 
product of the homogeneous spectral density and the correlation coefficient 
function: 


*pp ( ^l’ = 


%(“>P< r l’ r 2> 


(2.4-2) 
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When the correlation coefficient function is represented by Equation (2-3-4) 
in Section 2-3, the cross spectral density function of the excitation pres- 
sure will be: 


,-*■ -»■ 


^pp ^ r l ,r 2 = $ pp (“) t expC-AjK?) cos K?J [ exp(-A 2 Kr)) cos Kn] 


where 


(2.4-3) 


<j>pp (w) = the homogeneous spectral density of the excitation 

pressure 

B, = |x^ - x 2 1 = separation distance along x-axis 

t = |y-^ ~ y"2 1 = se P arat ic> n distance along y-axis 


(x^,y^) , ( x 2>y2^ = coordinates of and r 2 


A^,A^ = correlation decay constants 


is : 


2.4.1. 1 Four Edges Simply-Supported Rectangular Panels 

For four-side simply supported panels, the normal mode 


F. k = sin sin lEIZ 

d k Z b 


(2.4-4) 


Substitution of Equations (2.4-3) and (2.4-4) into Equation (2.4-1) gives the cross 
spectral density of the generalized force : 


Z Z 


jkmn 


. / \ r r , iirx-i muxo 

- <f> W J } e exp(-A]K?)cosK£ sin — sin— dx-^dx 2 

o o Z Z 


• f b , . kiryi niryo 

. j j e exp(-AnKq) cosKq sm sin dy.,dy 0 

n z u K 1 2 


* *pp<"> V, 

where I and I represent the two double integrals respectively, 
x y 


(2.4-5) 
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These integrals can be carried out by using the following identify: 


I 


x 


£ £ 

r r 3ir x i mirxo 

J J e expC-A^K £) cos sin — — sin dx^dx2 

oo & - 


£ x x 

= J j exp [-A3K (x^ - -x .^) ] cosK (x^ 
0 o > 


x^) sin 


~£ 


mirxo 

sin - dx£dx^ 


£ x 9 

t f ~ jrXi. 1077x2 

+ j J exp [-A K (X2 - x^) J cosK (X2 - x^) sin sin dx^dx2 

00 £ £ 


= Il+X 2 


( 2 . 4 - 6 ) 


where and I2 represent the two integrals respectively in the above 
expression. 

The identity of the above expression can best be seen from the following 
figure. If the variable and X2 are represented on rectangular coordinates ox-, 
and 0x2, then the region of integration of the left-hand side of Equation ( 2 . 4 - 6 ) 
is the square OABC. This square is divided by the diagonal OB into two isosceles 
right triangles. Inside triangle OAB, x-, > x 2 , while inside triangle OBC, 
x 2 > Xj. On the diagonal OB, xj = X2. It is easily seen that integral . 1 -^ repre- 
sents the integration over the .triangle OAB, and integral I 2 represents the 
integration over the triangle OBC. Therefore the sum of I± and I 2 is the inte- 
gration over the square OABC. 
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Carefully carrying out the integrations in eq. (2-4-6) yields 


v >V r\t <v 

T = = T + T + I + i 

x -2 xll xl2 xl3 xl4 


+ I + i + i + i 

x21 x22 x23 x24 


(2-4-7) 


Similar calculation gives 


I = — - = x + i + i +x 
y 2 yll yl2 yl3 y!4 


+ I + I + I „ + I ,, 
y21 y22 y23 y24 


(2-.4-7a) 


Substitution of eqs. (2.4-7) into eq. (2.4-6) yield the normalized cross 
spectral density of the generalized force as: 


where 




jkmn 

u jkmn - ^^^2 ^ “ ■ k xll ' J -xl2 


7 = ■■ -U-- — T +T -fT +7 +7 +X 

1 *t O -A- -I -I J- T o .o -*--1/ -Lrt-1 


xl3 x 14 x21 x22 


PP 


+ I „ + i +x +x + t „ + T + i 
x23 x24 yll yl2 yl3 y!4 y21 

+ I + I + I 
yll y23 y24 


(2-. 4-8) 


^ x11 S^ 2 + (mn+A l ) 2 J 

*[ I xlll ' X xll2 + 2 <- 1 > J+l <«W\> e ^ I^J 

X = i 

Xl2 S^ 2 + (nut-A 1 ) 2 J 

* [ I xl21 - ‘"'V I xl22 + 2< - 1)J+1 e ^ ^ 123 ] 

1x13 = Mm-J,) 2 ] [ Ix131 + 1x132 ‘ ^ Ixl33 ] 
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? w 


^“^Uv 2 ] ^ 111142 ■ 2a * 

1x21 8 [°1 2 + « w \) 2 ] 

X L^ll - «” + V ^212 + 2 <- 1 >”" 1 «” + I x213 ] 


x22 


x23 


x24 


S^ 2 + ( 3 *-\) 2 ] 

* [^21 - < J ” + V 2 x 222 + 2 <- 1 )'* 1 o * - V ^ 1 I x223 ] 

8^ 2 + (j!t-Ai) 2 J [ I)t231 + <J * ’ ? 'l ) ' 2C! 1 X x233] 

^ 2 7 o*^*] [W - <JI + V 2 x 242 - 2 “l W] (2 - 4 ' 


yii 


s[a 2 2 + (n*+A 2 ) 2 J 


X [lym - * V V 12 + 2t - 1)kfl + 4 2 > e ^ 113 ] 


.k+1 


712 8[q: 2 2 + ( nJ t-A 2 ) 2 ] 

X [l yl21 - <n* - A 2 ) I yl22 + (nn - I yl J 

JyU “ S^ 2 7 t Iyl31 + ( " ' " 2> W ‘ V33] 

^ “ 8 |a, 2 ♦ (!«4 2 ) 2 J [W ‘ (n ” + V V42 - 2a 2 S«s] 
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Ey21 s[a 2 2 + (k* + A 2 ) 2 ] 

x [w - < k ” + A 2 


n+1 ~^2 

x Vii - <k ” + v w + 2( - iy (k * + V e VuJ 


[y22 8[ct 2 2 + (kit - A 2 , 2 ] 

x (>21 ' (k ” ■ A 2 


r n+l -Q; 2 1 

x [V 21 - (kn - V V 22 + 2 <^> < k * ' V e I y 223 j 


y23 8 a 2 + (kK - A 2 ) 


>j [*y231 + (kit - A z ) I y232 -20^ I y233 ] 


8 a 


> ^7 i[w ' <k * + v w -> V43] (2 - 4 - 10) 


>11 = (-V^ 1 «L [(Ji)* * 2 A 2 + (-J)” + 2 \] sin ^ 

1 _ f_l)j +In l _ 1 ■ (~ 1 ) nrf ' 3 cos 2 \ 

^xll2 (j+m)it (j _ tn)it (nri'jjtt +■ 2A^ 

-1 + cos 2 A 

-J- 1 2 

(m-j)it + 2A^ 

a. sinA. + (jit + A^ cosA 1 -ql sinA^ + (jit - A) cos A. 

I xU3 = — \ \ + — J 1 2 (2.4-lla) 

+ (j<t + \y a L + (jit - AjT 

r 

>21 “ t' 1 ) 2 “l (j+m)it - 2?^ ’’ (j-m)it + 2^ sln 2A 1 

i - , -i + f-n-t-J t 1 - t- 1 ) 2 *” > 

xl22 (m+j)rt (m-j)Tt (j+m)Jf - 2A^ 

1 - cos 2A 1 

(j-tn)Jt + 2A 1 


17 



0^ sinA^ + (j« + A^) cosA^ 

r xl23 = 2 , . . ~2 

c* L + (j* + \) 


■a. 


sinA^ + (jit - Ap cosA^ 

a L 2 + (J* + \) 2 


x!31 


(-l)^ +m a. 


1 [( j + m) * - 


2 ( j -m) rt + 2 A^ 

j+m 


sin 2A^ 


-i + , i - f-iV"J , 1 ' <-*> c ° s 2 \ 

x!32 (nri"j)rt (m-j)it ( j+m) it - 2A^ 


1 - (-l)- 3 " 111 cos 2 A^ 
(j-m)Jt + 2 Aj^ 


(l) 3 ^ cos A^ + (j3t+ A^) sin aJ 

1 1 O O O <1 


-a 

e - a. 


x!33 


of + a«+ \) 

(-1) 3 ^ jct:^ cos A^ + (jit - 7^) sin A^J e ^ + 0!^ 


-a. 


a L 2 + (jit - Ap 2 


xl41 ^ ^ °1 (m+j)jt + 2A^ (m-j)jt + 2A^ S ^ n 2 \ 


-1 + f-l) 3+m . -1 + f-l') 3 ' m ,:L “ ( ~ 1 ^ nrl " :5 cos 2 \ 

^x!42 (j+m)tf (j-m)jt + 2A^ 


+ -1 + (-l) m 3 cos 2A^ 
(m-j)n + 2A^ 


(2.4-llb) 


(2.4-llc) 
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(-1) 2+1 a L cos A L + (j« + A^ Sin aJ 


-a. 


- a. 


x!43 


a L 2 + (j* + A x ) 2 


i+i r .1 ~ a i 

(-■1) J j cos A^ + (jn - A^) sin A^ e + 


a x + (j« - \y 


(2.4-lld) 


I x211 (_1)j a i [(j+iii)* + 2A 1 + (j-m)rt+ 2?^ sin 2A 1 

1 - (-1) 1 - <-l) m_J 1 - (-l) J+m cos 2Aj_ 

^x212 (m+j)rt (j +m )* + 2A^ 


-1 + (-L) j_m cos 2 A l 
* (j-m)rt + 2 A l 


sin A^ + (mJt + A ) cos A^ 

C x213 = „ 2 4 . / „ j. 7 .2 

+ (mrt + A^) 


sin A^ + (niJt - A^) cos A^ 

a ± 2 + (m* - A x ) 2 


(2. 4- lie) 


I x221 ( ' L ^ j a i [(nri-j)Jr - 2A^ + (m-j)* + 2 aJ Sin 2A 1 

nrf n 

_ 1 - f-l) j+m t -1 + f-l) j ~ m t 1 - ( - 1:> cos 2A 1 
x222 (j+m)rtn (j-m)Jt (nri'j)rt - 2A^ 

1 - (-l) m " j cos 2 A l 

(m-j)« + 2 Aj^ 
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sin A^ + (mit + A^) cos A ] 

E x223 2 7 9 

0!^ + (mJt + A ) 


-a^ sin A^ + (mJt - Ap cos A^ 
« L 2 + (m* - A l ) 2 


I = a 

x231 1 


: 1 + 1 1 . 

1 (ro+j) 11 - 2A^ (m-j)x+ 2A^ sin 


2 A. 


- -1 + C-l) j+ ° , 1 - <-lv»"° , 

f *-J \ ~ ✓ . x ^ *»* 


_ f \ j “ m 1 - (-l) 1 ”^ cos 2 A. 


x232 (j+m)n 


(j-m)it - 2 A 


1 - (-l) m " J cos 2A 1 
(m- j ) it + 2A 


x233 


HI+ i F ~i ~ 

(-1) J^- cos A^ +■ (nut + A^) sin A^ I e ~ a i 


Ct^ + (mtt .+ A-)' < 


m +i n ~j ~ o** 

cos A^ + (mrt - Ap sin A^J e + 


Q^ 2 + (mJt - A^) 2 


W = “ 


* r 1 

1 |^(j+m)n 


+ 


+ 2 \ (j-m)« + 2 A l 


sin 2 A, 




- 1 \ m- 3 1 ” C-l) j m cos 2A 


x242 (m+j)Jt (m-j)jt (j+m)jr + 27^ 


-1 + (-l) J " m cos 2A 1 
(j-m)x + 2 A 


(2.4-llf) 


(2.4-llg) 
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*x243 “ 


+ 


I ylll “ 


I yll2 


+ 


X yll3 


+ 


X yl21 


x yl22 ° 


+ 


(„l) in+X cos X]_ + (m-rr + A-^) sin A., ] e X - a-j 

- 2 

a-^ + (urn + A^) 

nH~l f" "1 — 

(- 1 ) (_ a l cos a i + (nnr ~ A-^)sin A^ J e + 

2 2 

0 £^ + (m?r - A-^) 


(2.4-llh) 


.n+k 

(-1) a. 


r — - — + — ~ — isi 

L (n+k)fr + 2Ao Cn-kW + 2A«-J 


(n+k) it + 2A 2 (n-k)ir + 2A 2 * 
k-n 


sin 2\, 


1 - C-l)^ + 1 - + 1 - M) n+k c°s 2X 2 

(k+n)ir (k-n)ir (n+k)7r + 2 A 2 


-1 + (-l) n k cos 2A 2 
(n-k)ir + 2 A 2 

a 2 sin A 2 + (k?r + A 2 ) cos A^ 
a 2 2 + (kir + A 2 ) 2 

~ a n sin A 9 + (kir - A,) cos A~ 

— f ~ ^ (2. 4-12a) 

+ (ku - A 2 ) 

(-W a 2 [ (k+n)» - 2,1 2 + (k-n)it + 2 ~ 2 J sln 2X 2 

1 - (-l) n+k + -1 + (-l) n - fc + 1- (-l) k+n cos 2\y 

(n+k) u (n-k)ir (k+n) it - 2X 2 

1 - (-l) k_n cos 2A 2 
(k-n)ir + 2A 2 
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CC 2 sin + (k + A^) cos A^ 

r yl23 = 2 7~2 

^ a 2 + (kJt+ A 2 ) 


-a 2 sin A 2 + (kjt - A 2 ) cos A 2 
a 2 2 + (kit - A 2 ) 2 


(2. 4- 12b) 


Vsi ~ (-i) k+n 


' i_ 

(k+n)jr- 


2A 2 (k-n)* + 2A 2 


sin 2A„ 


,k+n 


,n+k 


l -1 + (-l)"' 11 | 1 ~ (-l) k " n , 1 ~ ( ~ 1) C ° S 2A 2 


y!32 (k+n)jt (k-n)jT 


(n+k) jt - 2 A 


1 - (-l) n " k cos 2A 2 
(n-k>t + 


k+1 T 1 -Ct 2 

(-1) J^-or 2 cos A 2 + (lot + A 2 ) sin A 2 J e -a. 


yl33 


a 2 2 + (fat + A 2 ) 2 


k+l T T *"®2 

(-1) a cos A + (kit - A ) sin A e -a, 

L ^ ^ i 


ct 2 2 + (tot - A 2 ) 2 


(2.4-12c) 


X yl41 “ (~ 1 > n+k «2 


" 1 . 

(n+k)rt 

_ 


+ 2A 2 (n-k)jt + 2A 2 


sin 2A„ 


i s j- Tr 

-1 + (-l) k+n -1 + (-l) k " n 1 " (-1) COS 2A 2 

X yl42 (kt*n)Jt + (k-n)rt + (n+k)rt + 2A2 


-1 + (-l) n " k cos 2A 2 
(n-k)« +' 2A 2 


22 



y!43 


(-l) k+X cos + (kit + A^) sin A J e - Ct. 


a 2 1 2 + (kit + A 2 ) 2 





cos A^ +■ (kit - A^) sin A 2 
a 2 2 + (kit - A 2 ) 2 



cl 


I y211 a 2 [(Id-n)Jt + 2A 2 + (k-n)rt + 2 aJ sin 2 \ 

i - (-n n+k 1 - (-n n - k 1 ' <~ 1)k+n cos 2A 2 

y212 (n+k)it (n-k)it (k+n)it + 2A 2 


-1 + (-l) k n cos 2A 2 
(k-n)it + 2A 2 


X y213 


sin A 2 + (nit + A ) cos A 2 
« 2 2 + (n< + A 2 ) 2 


-a 2 sin A 2 + (n« - A 2 ) cos A 2 

4- ■ i — , 

a 2 2 + (nit - A 2 ) 2 

X y221 = (_1) a 2 [(n+k)Jt - 2A 2 + (n-k) it + 2 aJ sin 2 \ 

i . c-i) k+n -i -f (-i) k " n 1 ' ( _1 ) n+k cos 2 \ 

X y222 (k+n)it + (k-n)it + (n+k)it - 2A 2 

1 - (-l) n ~ k cos 2 A 

4- ■ — -• 

(n-k) it + 2A 2 


(2 .4-12d) 


(2.4-12e) 
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a 2 sin A 2 + (n« + A 2 ) cos 

I y 223 = T 2 o 

2 + (n?c + A 2 ) 


-a 2 sin A 2 + (nit - A^) cos A 2 


a 2 2 + (nit - a 2 ) 2 


X y231 (-1) a 2 [(n+k)it - 2A 2 + (n-k)fl + 2 aJ sin 2A 2 

n+if 

-1 + (-l) k+n 1 - (-l) k_n 1 ‘ ^ _1) cos 2A 2 

^y232 (k+n)it (k-n)it + (n+k)it - 2A 2 

1 - (-l) n " k cos 2A 2 
(n-k)it + 2A 2 


y233 


(_l) n ^ -CX 2 cos A 2 + (nit + A 2 ) sin A 2 J 


-a. 


- OL 


2 2 
a 2 + ( nlt + A 2 ) 


n+1 P 1 ~ a 2 

(-1) a 2 cos A 2 + (nit - A ) sin A 2 e + a, 

i '■ ' 


a 2 2 + (nlt - A 2 ) 2 


^241 a 2 [(k+n) it + 2A 2 + (k-n)it + 2A 2 J Sin 2 \ 

-l + (-i) n+k -l + (-n n “ k 1 _ (_i) cos 2 \ 

^242 (n+k)it + (n-k) it + (k+n)it + 2A 2 


>1 + (-l) k ~ n cos 2A 2 
(k-n)it + 2A 2 


(2.4-12f) 


(2. 4-12g) 
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I y243 + 


n+1 T 1 -a 2 

(-1) I " a 2 cos ^2 + ( nir + ^2^ sin *2 J e ~ a t 

2 “ ~2 
+ (nir + 


(-1) £ &2 cos ^2 ( n1T _ ^2^ s ^ n 


].■“* + « 2 


a 2 4- (mr - A 2 ) 2 


(2 . 4-12h) 


Note: When j=m or k=n, the indeterminate terms (5.) in the above expressions 

for the I should be set equal to zero- ^ 

= A-jKl 

a 2 ~ ^Kb 
X 1 = k* 
x 2 = Kb 

K = - (2.4-13) 

c 


mode is: 


2. 4. 1.2 Four Edges Clamped Rectangular Panels 

For the four edge clamped rectangular panel, the normal 


V - x j fe) \ (Y > 


where 


X^(x) = cosh x-^ - cos x^ 

- 0.9825 (sinh x-^ - sin x^) 
Y (y) s= cosh y^ - cos y^ 


- 0.9825 (sinh y-^ - sin y^) 

X. (x) = cosh x. - cos x- - (sinh - sin x.) 

3 J J 3 3 

3 = 2, 3 ... 
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Vy> 


cosh y fc - cos y k - (sinh y k - sin y k ) 


k = 2, 3 . . . 


- 

1.5056 ttx 

X 1 “ 

Z 

*1 = 

1.5056 ty 
b 

- 

(.1 + 1/2) rrx 

X • 
1 

Z 

*k = 

(k + 1/2) rry 

b 


(2.4-15) 


For this case the carrying out of the integration in Equation (2.4-1) is too 
complicated and is not practical. Therefore the calculation for the has 

to be done by numerical integration. Substitution of Equation (2.4-15) into 
Equation (2.4-1) and representing the double integration by double summation 
yields 


1 j kmn ^ 
$ pp (w) 


Z Z 

J J exp^AjKg) cos(Kg)x. (x x ) ^ (x 2 ) dx 1 dx 2 
oo 

b b 

. / / exp(-A 2 Kq) cos(Kp)Y k (y x ) Y n (y 2 ) dy 1 dy 2 
o o 


= II exp (-AqK£) cos (K^)Xj (x p ) ^ (x q ) Ax]_ Ax 2 

p q 

• I I exp(-A 2 Kq) cos(Kr})Y k (y p .) Y n (y ,) Ay 1 Ay 2 

p' q' 

When the length Z and width b are divided into s and s’ equal divisions re- 
spectively, simple calculation gives 


Ax^ = 

Ax. 

Ay x - 

Ay ; 


Z 

s 

b 
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1/21 

1 / 2 ) 

= 2, 3 

= 2 , 3 (2.4-17) 

(2.4-16) yields the normalized 
as 




(2.4-18) 


• “*[“ Cp - ] x j \ 


' exp [~ |p' - q' | ]c°s[|t (p' - q') ] Y k Y n 


where the X and Y functions are given by Equation (2.4-15) with the x and y 
arguments given by Equation (2.4-17). 

The numerical method described above can be applied to any correlation 
function and any mode shape . 

2.4.2 Gross Spectral Density of Generalized Force, Exponentially 
Decaying Correlation 

When the correlation coefficient is represented by Equation (2.3-3) 
in Section 2.3, the cross spectral density function of the excitation pressure becomes 


4> pp (r 1 ,r 2 ,m) = $ pp (m) exp(-AqKg) exp (-A 2 Kp) 


(2.4-19) 


Substitution of Equations (2.4-19) and (2.4-4) into Equation (2.4-1) yields 

£ £ 

// 

o 0 


I jkmn = / / ex P (“ A 1 K S) sln sin dXqdx 

oo £ £ " 


b b 


kiry 


t r n wy 9 

/ J exp(-A 2 KT]) sin — — sin — g— =• dyq dy 2 (2.4-20) 


o o 


The above integration can be carried out by the following identity: 


£ £ 


/ / exp (~A]K£) sin sin dx dx, 


o o 


12 


^ i j-n-Xq Ttnrx 2 

= / dxq / exp [AqK (xq - x 2 )]sin — - — sin — - — 


<1 dx 2 
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£ X 2 jrrx^ m7rX 2 

+ / dx 2 / exp t-A 2 K(x 2 - x^Jsin — g — sin — ^ — dx^_ 

o o 


= *1 + *2 


(2.4-21) 


Carrying out the integrals 1^ and I 2 yields 


I j m 


£ £ 


H A/ A/ 

= — 2 ’ I J exp (“A K£) sin 

P -** 


0 0 


17tx 2 mirx 2 

— - — sin dx, dx 

£ £ l 2 


r i , i i 

2 L cl ’ 1 + ( j it) 2 a-, 2 + (m*) 2 j 


2 + [ (-l) j ' +1 + (- 1 ) Tttfl ] exp (-«!) 

+ (itiXuitt) 2 2 ir 2 zl for j = m 

“l + (;iir) J I «2. + (mu') J 

(2.4-22) 


jL £ £ jirx^ nnTX2 

= ~2 / / exp (-A-jKTI ) sin sin — ^ — dx^dx 2 

& _ . . 

O 0 


(jir) (rait) 


2 + [ (-1) j+1 + (-l)”* 1 ] exp (-a 1 ) 
[ c ^ 2 + (j-rr) 2 ] [ a 2 2 + (but) 2 ] 


j - r (-p m J -i t ( 

v-l + (jir) 2 L 2 2 

_r (-p j 

( m -iT) 2 L 2 


m+j 


2 (m+j) 

j+tn 




a 2 + 


- 1 + (“D 


2 (j+m) 




for j 4 m 


(2.4-23) 


Similarly the double integration on and y 2 in Equation (24-20) can be carried 
out to give 
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L kn 


b b 


1 r kfry 1 niry 2 

2 — / / exp(-A 2 Kn) sin ± sin _ dy 1 dy 2 


o o 


a i r i 

2 2 4 

La. + 


(kir) (*2 + (mi ) 2 


] 


kn 


2 + [ (“D + (- 1 ) exp (- 02 ) 

(kir) (nir) — 2 ^ for k = n 


^a 2 2 + (kir) 2 J£a 


, 2 + (nir) 2 ] 


(2.4-24) 


b b 


Iff , „ . niry 2 

— j J ) exp(-A 2 Kri) sm — — sm — - - — dy^dy 2 


o 0 


(kir) (nir) 


2 + [(-D k+1 + (-l) n+1 ] exp(-a 9 ) 


|^a 2 2 + (kir) 2 J^a 2 2 + (nir) 2 ] 


2 , .2 

a 2 + (kir) 


[ , n-k n+k 

(-D “1 4. (-D -1 I 

2 (n-k) 2 (n+k) J 


n 


a 2 2 + (hit ) 2 


[ 


C-l) k “ n -l (-l) k+n -i 
2 (k- n) 2 (k+n) 


] 


for k ^ n 

(2.4-25) 


Substitution of Equations (2.4-22), (2.4-23), (2.4-24), and (2.4-25) into Equation 
(2.4-20) yields the cross spectral density of the generalized force of the ex- 
citation as : 


I jkmn 


(b£) 2 

$ 

PP 


I jm I kn 

(bf) 2 

*PP 

(w) 

T T^* 
jm kn 

(b£) 2 

$ PP 

(to) 

a 

T T 

jm km 


for j - m, k = n 

for j = m, k ^ n 
for j ^ m, k = n 
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= (b£) 2 4> pp (oj) Ij m 1^ for j m, k ^ n (2.4-26) 

2.4.3 Joint Acceptance and Normalized Cross Spectral Density of 
Generalized Force 

The normalized cross spectral density of the generalized force 
of the excitation is defined as 


where 


■'■jkmn 


■^jlcmn 
S 2 * pp (w) 


S = b£ = area of the panel 


(2.4-27) 


The relation between the joint acceptance squared and the normalized cross 
spectral density of the generalized force is given by [see Reference 3, Chapter 
VI, Section 6.2, Page 24, Equation (6-25) and Reference 2, Chapter 7, Section 
7.2, Page 212]. 


J7, 

jkmn 

I..,-, COS “ (t,-, ~ t ) 

jkmn v jk mn 


(2.4- 

II 

d 

phase difference of the jk and mn 

modes 


kmn 

[l- («A> jk ) 2 ] [ 1 - (“/“mn) 2 ] + 

2 

4t-vt m „ at /(u., o ) 
^jknnn jk mn 7 

^jkmn 

A?. + B 2 

jkmn jkmn 



^ j kmn 

2 |[l “ (“ /£0 jk )2 ] ? mn“ /a) mn " [ X " 

(m/a)) 2 

mn 

] V'V } 


cos u(T jk - x mn ) 


^•j kmn 
^jlcmn 


(2.4-29) 
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2.5 The Displacement Response 

The displacement response spectral density at any point r(x,y) is given 
by [see Reference 3, Section VI, Equation 6-24, Page 24]. 


^ww ( r »u) 


where 

$ (r.co) 
ww v ’ 




l 

j,k,m,n 




|h* |J.i , 

1 mn 1 j kmn 


(2.5-1) 


displacement spectral density in (inch) 

rad . / sec . 


S’ = b’*’ 


area of panel subjected to excitation in (inch) 


b' = width of panel subjected to excitation in inches 


= length of panel subjected to excitation in inches 


F jk’ F - 


mn 


normal modes 


$ (w) 

PP V ' 


spectral density of the spatial homogeneous excita- 
tion pressure in (psi) 2/radian per second 


The magnitudes of the complex frequency response functions and their complex 
conjugates are given by 


where 


V 1 

II 

J* 

[‘“V ■ 

2,2 

w ) + 

( 2 ? jk u jk t °) 2 

j - 1/2 

H* 1 
mn 

- “nm 1 


2.2 , 

10 ) + 

2 

( 2 L.« «) 
ran mn 

j - 1/2 

(0 

= frequency in 

radians 

per second 



- 2TTf 


(2.5-2) 


f = frequency in Hertz 
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The modal masses are 



«mn 



dr 


M?-b 

4 


with 


(2.5-3) 


M = ph + p*h' = mass distribution per unit area 

The joint acceptance squared for different mode combinations are given by 
Equation (2.4-28), in Section 2.4. 

2 

The displacement response spectral density in inch /Hertz is given by 


S w< r -« 


2u$ (r.to) 


(2.5-4) 


o 

The mean-square displacement response in inch is given by 

w2(?) = / S w £,f)df 

The root-mean-square displacement in inches is given by 
w(r) = tw 2 (r)] 


(2.5-5) 


(2.5-6) 


2.6 The Excitation Pressure Data 


For derivation of the formulas, see Reference 4, Section 6.1. The excita- 
tion spectral density in (psi) 2 is given by 

radian/ sec. 

Spp(f) - 170.576 

jo 

V (<o) = W [1 °] (2.6-1) 

The excitation spectral density in decibels/Hertz is given by 

Sp p (f) = S 3rd (f) " 10 los 10 (0* 23157f ) (2.6-2) 
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where 


Sg^Cf) = one- third octave excitation pressure level 

in decibels (input) 

O 

The excitation spectral density in (psi) /Hertz is given by 


S' (£) = 2tt* (to) 

PP PP 


(2.6-3) 


The relations between the overall mean-square pressure and the overall pressure 
level of the excitation are: 


p 2 = 8.75526 x 10“ 18 [10] La/1 ° 

3 . 


(2.6-4) 


and 


L a = 170.576 + 10 log 1Q (p a ) 


(2.6-5) 


where 


2 2 
p = mean-square pressure in (psi) 


L a = excitation overall pressure level in decibels 
referenced 0.00002 Newton/meter 2 

2.7 The Acceleration Response 

The acceleration response spectral density at any point r(x,y) in 
/ inch \ 2 1 . 

(^7 ) 1S Siven by 


A + 

*ww( r »") = “ *sw( r »“) 


ww' 


(2.7-1) 


where the displacement response (^(rjW) is given by equation (2.5-1) in 
Section 2.5. 


The acceleration response spectral density in g Z /Hertz is given by [see 
Reference 4, Section 6.1, Page 18], 

S-(r,f) - 4.215093 x 10" 5 *-(£,») 


(2.7-2a) 


The acceleration spectral density in decibels referenced u g" 


S wg (r ’ f) = 10 iog 10 [S w (r ’ f)3 


(2.7-2b) 
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The mean-square acceleration response in "g 2u is given by 


G 2 £) = / S^(r,f)df (2.7-3) 

The root-.mean square acceleration in "g" is given by 

G(r) = [G 2 (r)] 1/2 (2.7-4) 


2.8 The Stress Response 


For derivation of the formulas, see Reference 3, Section 6.1. The stress 
response spectral density in (psi)2/radian per sec. is given by 


l> ao (r,m) = y (r)$ w (r,a>) 


( 2 . 8 - 1 ) 


where the displacement response spectral density ^^(rjW) is given by 
Equation (2.5-1). The factor to convert the displacement response into 
stress response is given by 


2 ■* 
Y (r) 


(Eh,) 


4(l-v 2 ) 


2 2 
Q x + Q y 

Q 2 w 


( 2 . 8 - 2 ) 


h, = h + hr 


largest height of stiffeners at r (See Figure 1) 


^x = 


m,n 
=1,3 . 


TT 

mn 


[(f ) 2 + <^f] 


mrx . nmy 
sin — r- sin r~ 

X, D 


D x(f) + 2H (ff) + MD 

(2.8-3) 


Q y " 


I 

m,n 
= 1,3, 


»4[(?) 2 + v(=) 2 ] 


. mirx . nry 
sm — - — sin - 

x. b 


(?) - + M?) 


X 


(2.8-4) 
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0 - 2 
tj m,n 

' X j 3 J at • 


mmx 


sin. 


sin 


liny 


mn 


[MS) +2H (^) + M?) ] 

The stress response spectral density in (psi)^/Hertz is given by 


S CTa (r,f) = 2ir« (r,u) 


n2 


The mean-square stress response in (psi) is given by 


9 

u (r) 


/ Sera (r ,f )df 


2 

The root-mean-square stress response in pounds/inch is given by 


U- 


( 2.8 


( 2.8 


a(x) 


[c 2 (v )] 1/2 


( 2.8 


2 . 9 Average Responses Over the Whole Structure 


The average response over the whole structure is obtained by integrating 
and dividing the result with the area of the structure. Detail derivation of 
the formulas is given in Reference 5, Section 2.5. The final results are 
summarized, as follows: The average displacement spectral density in inch”/ rad. 

sec. is 


$ (m) = 

T.TT.T V * 


WW 


S 4> PP ^ J,. I 3C I y f H jlc I J jkjk 


1 >k 


y' jk 1 


(2.9 


where the expressions for I x and I for different boundary conditions are as 
follows : 

2.9-1 Four Edges Simply-Supported Panels 


I 

x 



(2.9 


2.9.2 Four Edges Clamped Panels 


I (1) = 

x 


I (1) = 


1,5056m ^ (l*5056)ir- 1 + sinh 2(1. 5056m) 


-sinh 2 (1.5056m) - 2 exp (-1.5056m) 


(2.9 


8-5) 

- 6 ) 

7) 

- 8 ) 

/ 

- 1 ) 

- 2 ) 

-3a) 
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I x ~ I x^ ~ — [ (j + ^)ir - 1 + \ sinh 2 (j + -|)Tr 

(j + p ir L 

-sinh 2 (j + -j)tt “ 2(-l) 2+1 exp [-(j + -|)?r]j 
j = 2, 3, ... (2.9-3b) 

Iy = I (k) = f + .jp 1 * ~ 1 + ~2 sinh 2 (k + ^r)ir 

(k + ^ir 

-sinh 2 (k + ^) - 2(-l) k+1 exp [-(k + -|) 

k = 2, 3, ... (2. 9-3c) 

2.9-3 Two opposite edges simply-supported while other two clamped panels 

X x = \ (2.9-4) 

I^(k) is given by Equations (2.9-3a) and (2.9-3c) 


2.10 One-nth Octave Frequency Increment and Number of Data Points 

In the calculation of the response spectrum as a function of frequency, 
uniform frequency increment is not convenient, because when the spectrum is 
plotted with the frequency in logarithmic scale, the points will be too close 
in high frequency and too separated in low frequency. In order to obtain good 
plots and save computer time, it is the most convenient to use cne~nth octave 
frequency increment in the calculation. 

The interval from frequency to f -2 will be one-nth octave if 


f \ n 

4 ) ■ 


or 


= 2 


1/n f 


(2.10-la) 


(2.10-lb) 


The bandwidth of the one-nth octave band expressed in terms of the lower limit 

(2.10-lc) 


f-^ is: 


Af 


f 1 = (2 1/n - 1) f x = D^f 


hi 
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where 


1 

D' = 2 n -l (2.10-ld) 

n 


The geometric mean frequency of f-^ and f ^ is 

f = (f x f 2 ) 1/2 (2.10-2) 

Solving equations (2. 10- la) and (2. 10- lb) gives- the lower and the upper 
limits of the one-nth octave band in terms of the geometric mean frequency: 

- 1 

f x = (2) 2n f 

JL 

f 2 = (2) 2n f (2.10-3b) 

Therefore the band-width of the one n-th octave band is 


Af = f 0 - f-. = D f 
2 i n 


(2 . l0-3c) 


where the one n-th octave band width constant is given by 



- 1 
"2n 


(2.i0-3d) 


The values of and D n for n from 1 through 50 are as follows. Note that 
when n equals to 1, 2, and 3, the values of and D n are the corresponding one- 
octave, one-half-octave and one-third-octave band width constants. 
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Value of 


Value of n 



Value of 

l/2n -i/2n 
D = 2 -2 


1 

1 . 

0.707107 

2 

0.41421 

0.348311 

3 

0.25992 

0.231563 

4 

0.18920 

0.173504 

5 

0.14869 

0.138740 

6 

0.12246 

0.115588 

7 

0. 10409 

0.0990615 

8 

0.09050 

0.0866705 

9 

0.08006 

0.0770354 

10 

0.07177 

0.0693286 

11 

0.06504 

0.0630238 

12 

0.05946 

0.0577703 

13 

0.05476 

0.0533253 

14 

0.05075 

0.0495156 

15 

0.04729 

0.0462139 

16 

0.04427 

0.0433251 

17 

0.04161 

0.0407762 

18 

0.03925 

0.0285106 

19 

0.03715 

0.0364835 

20 

0.03526 

0.0346591 

21 

0.03355 

0.0330085 

22 

0.03200 

0.0315080 

23 

0.03059 

0.0301380 

24 

' 0.02930 

0.0288821 

25 

0.02811 

0.0277268 

26 

0.02701 

0.0266603 

27 

0.02600 

0.0256728 

28 

0.02506 

0.0247559 

29 

0.02419 

0.0239022 

30 

0.02337 

0.0231054 

31 

0.02261 

0.0223601 

32 

0.02189 

0.0216613 

33 

0.02122 

0.0210048 

34 

0.02059 

0.0203870 

35 

0.02000 

0.0198045 

36 

0.01944 

0.0192544 

37 

0.01891 

0.0187340 

38 

0.01840 

0.0182410 

39 

0.01793 

0.0177732 

40 

0.01748 

0.0173289 

41 

0.01705 

0.0169062 

42 

0.01664 

0.0165037 

43 

0.01625 

0.0161199 

44 

0.01587 

0.0157535 

45 

0.01552 

0.0154034 

46 

0. 01518 

0.0150686 

47 

0.01485 

0.0147479 

48 

0.01454 

0.0144407 

49 

0.01424 

0.0141460 

50 

0.01395 

0.0138631 
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When the one nth octave bandwidth is used for frequency increments, the expression 
for the number of data points can be derived as follows. 

By definition of the one-nth octave band, the number of band between frequencies 
f and 2f is n. Hence for the geometric series 


f 


o 


f 


1 


f 


2 


f . 
1 




(2.10-4a) 


between each consecutive pair of frequencies, there are n bands. Therefore the 
total number of data points is 


i 


But 


2 J 


or 


3 


nj + 1 


l°g -f 

2 £ 

log 2 


Substituting into Equation (2.10-4b) gives 


(2.10-4b) 


f . 

log y~ 

i - n (2.10-4c) 

log 2 

When the result of Equation (2.10-4c) is not an integer, it should be increased to 
the next integer. Equation (2.10-4c) is very convenient for hand calculation to 
determine the number of data points. For example 


n - 33 

fo = 5 

- 5000 


40 - 



Substituting into Equation (2.10-4c) gives the number of data points 


i = 33 


, 5000 

lOg — r“ 


log 2 


+ 1 - 331 


41 



III. ANALYSES OF THE RESULTS 


3 . 1 Effect of Boundary Conditions 

For a typical run (Run No. 8010) of the Program RANDOM with input data 
described in Section 5.7, the frequency of the fundamental mode, the displace- 
ment, the stress and the acceleration spectral densities at fundamental mode, 
the root-mean— square displacement, stress, and acceleration at the center of 
the panel for the three boundary conditions are as folloi^s:' 


Boundary 

Condition 

Fundamental 

Mode 

Frequency in 
Hertz 

Displacement 
Spectral 
Density at 
Fundamental 
Mode in in. / 
Hertz 

Stress Spectral 
Density at Funda- 
mental Mode in 
(psi) ^/Hertz 

Acceleration 
Spectral 
Density at 
Fundamental 
Mode in g2/ 
Hertz 

SSSS 

148.0 

0. 48339292E-05 

7304.0415 

24.259805 

CCCC 

215.2 

0.37594843E-05 

5680.5610 

84.278543 

SCSC 

206.4 

0 . 18984536E-05 

2868.5533 

36.038467 

Here 





SSSS 

= Four edges 

simply-Supported 



CCCC 

= Four edges 

clamped 

— 

— 

SCSC 

= Two opposite edges simply-supported while other 

two clamped 

Boundary 

RMS 

RMS ‘ 

RMS 

RMS 

Condition 

Displacement 

Stress in 

Acceleration 

Excitation 


in inch 

psi 

in g 

psi 

SSSS 

0.010248574 

398.37757 

23.549822 

0.22680647 

CCCC 

0.011197214 

435.25270 

49.516272 

0.22680647 

SCSC 

0.0078521891 

305.22652 

34.075105 

0.22680647 


It is interesting to note the following results from the effect of the 
boundary conditions: 

(1) The more rigid the boundary condition, the higher the natural frequency. 
Hence, the four edges clamped panel has the highest frequency, the SCSC panel next 
and the four edges simply-supported panel has the lowest frequency. 

(2) The more rigid the boundary condition, the less the displacement and 
the stress spectral densities at resonance. The panel with four edges clamped has 
the least displacement and stress spectral densities at the fundamental mode, the 
SCSC panel next, and the panel with four edges simply- supported has the largest. 

This conforms to common sense reasoning. 
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(3) The four edges clamped panel has the largest acceleration spectral 
density at the fundamental mode and the largest root-mean-square acceleration, 
the SCSC panel has the second largest, while the four edges simply-supported 
panel has the lowest. This is because that acceleration is proportional to the 
square of the frequency, and acceleration spectral density to the fourth power 
of the frequency. 

(4) The root-mean-square displacements and stresses for all three cases 
are not much different. The four edges clamped panel is 10% higher than the four 
edges simply-supported panel. The SCSC panel is 20% lower than the four -edges simply- 
supported panel. 

The plots of the response spectral densities are shown in Figures 2 through 

10 . . 


3.2 Contribution of Each Mode 


The contribution of each mode to the responses can be estimated from the 
response spectral density at each mode and the half-power bandwidth at each fre- 
quency. The contributions of each mode to the mean-square displacement and 
acceleration, are given respectively by 


w i = 

^i G ww 

(3.2-la) 

G 2 = 

B.S — 


X 

x ww 


where o 

a 

w. = 

X 

contribution of the ith mode to the mean-square 


displacement 


G i " 

contribution of the ith mode to the mean-square 
acceleration 


S ww = 

displacement spectral density at the ith mode 


S ww 

acceleration spectral density at the ith mode 


The half-power bandwidth of the ith mode is given by 


B. 

2S±fi 


i 

(3.2-2) 

where 

4 -= 

X 

damping ratio at the ith mode 


f i = 

frequency of the ith mode 


Here in the typical 

test run described in Sections 3.1 and 5.7, the damping 

ratio for all modes 

is 


4 

0.04 
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The contributions of the fundamental mode to the mean-square displacement and 
acceleration are calculated as follows: 


Boundary 

Condition 

Bandwidth At 
Fundamental 
Mode in 
Hertz 

Fundamental 

Mode 

Contribution ' 
to M.S. Dis- 
placement in 
inch2 

Fundamental 

Mode 

Contribution 
to M. S. 

Acceleration 

2 

xn 

% of M.S. 
Displacement 
Contributed 
by Funda- 
mental Mode 

% of M.S. 
Acceleration 
Contributed 
by Funda- 
mental Mode 

SSSS 

11.85 

57. 2 x 10 -6 

288.5 

54.5 

52 

CCCC 

17.22 

64. 7 x 10" 6 

1453 

51.8 

59 

SCSC 

16.55 

31. 41 x 1(T 6 

597 

51.0 

49 


From the above results, it is seen that the contribution of the fundamental 
mode to the mean-square response amounts to approximately 50% of the total mean- 
square value for both the displacement and the acceleration. 

It is to be noted that when the mean-square value contributed by the 
fundamental mode is used to estimate the root-mean— square response, the result 
will' amount to 


(0.50) 


1/2 


0.706 


70.6% 


of the actual root-mean-square response for both the displacement and the accel- 
eration. The detail values for each boundary condition are tabulated as follows: 



RMS 

RMS 




Displacement 

Acceleration 

Estimate 

Estimate 

Boundary 

Estimated ( 

Estimated 

Amounts 

Amounts to 

Condition 

from 

from 

to % of 

% of Total 


Fundamental 

Fundamental 

Total rms 

rms 


Mode in inch 

Mode in g 

Displacement 

Acceleratii 


SSSS 

0.00757 

16.95 

73.8 

72 

CCCC 

0.00804 

38.1 

71.8 

76.8 

SCSC 

0.0056 

24.4 

71.3 

70 


It is seen that the estimation of the root-mean-square response by the principal mode 
yields a result approximately 28% lower both in displacement and acceleration. 

The contributions of the (3,1) mode to the mean-square displacement and 
acceleration are estimated as follows: 
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Boundary 

Condition 


(3,1) Mode 
Frequency 
in Hertz 


Displacement 
Spectral Density 
At (3,1) Mode 
in in^/Hz 


Acceleration 
Spectral Density 
at (3 , 1) Mode 
in g2/Hz 


571.6 

764.7 
548.4 


9.848 x 10 -10 
5.254 x lCrJ-jj 
3.528 x lCf 10 


1.099 

1.879 

3.335 


Boundary 

Condition 


Bandwidth 
at (3.1) 
Mode in 
Hertz 


(3,1) Mode 
Contribution 
to M.S. Dis- 
placement in 

inch^ 


(3,1) Mode 
Contribution 
to M.S. 
Acceleration 


% of M.S 
Displacement 
Contributed 
by (3,1) • 
Mode 


45 x 10 9 
32.3 x 10" 9 
156 x 10“ 9 


0.0428 

0.0289 

0.252 


9 

4.67 

12.6 


It is seen that the higher mode contributes more to the acceleration than to 
the displacement response. 


3.3 Contribution of the Cross-Terms 

In the expression for the displacement spectral density, the summation 
may be separated into two groups of terms - the main terms and the cross terms 
as follows: 




= s' * (w) y 

pp l 

j ,k,m,n 

j=m,k=n 

+ S, *pp(«) l 

j ,k,m,n 
j t m 

k ^ n 


F jk< r > F mn< r )l H 3 k l|H nn |j Jkml (») 


V r) F »n( r )l H 3 kll H „nl J 3k» n (“) 


(3.3-1) 


where 


$ (r,io) 
ww 


= displacement spectral density at any point r as function 
of frequency to 


area of structure subjected to excitation 
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V Cco) 


F j k ’ F mn 


J 2 . 

jkmn 


Homogeneous spectral density of the excitation pressure 
normal modes 


magnitudes of the frequency response function 

joint acceptance of the jk and mn modes combination as 
function of frequency 


It is generally true that the cross terms contribute only a small portion 
of the total displacement. This is because | H .. | |H | is much smaller than 
• 1 2 ' J 11111 

IH.jJ . As to how small the summation of the cross-terms is and how it affects 
the spectrum of the response, there is limited information available in the 
literature. 


Program RSRPC3 is a modification of Program RSRPC1 designed to study the 
contribution of the cross-term3 and their effect on the response spectrum. 

Typical results of this program are given in the following. The input data 
for these results are the same as those of Section 5.7 except the radius of 
curvature is 70 inches instead of 100 inches. The fundamental mode frequency 
is 216.92 Hertz. Data in item 1 are the results of summation of all the terms 
while data in item 2 are the results of the summation of the main terms only. 

The range of the index j,k,m,n is from 1, 2,.,, to 5. That is 5^ = 625 terms are 
summed for each frequency to obtain the spectral density. The plots of the dis- 
placement and acceleration responses for four locations are shown in Figures 
11 to 34. It is seen that the cross- terms contribute not very much to the 
response but affect the response spectrum to some degree. It is to be noticed 
that the cross-terms sums up to positive values in some locations and negative 
in other locations. 


a. 

Location: Center of Panel 

1 

X “ 2" 

23.75 

inches 



b 

y = i - 

29.187 

inches 


Run No. 4220 Frequency 

Increment 

n =* 33 



Item 


Displacement 
Spectral 
Density at 
Fundamental 
Mode in inch 2 /' 
Hertz 


Stress Spectral 
Density at Funda- 
mental Mode in 
(psi) 2 /Hertz 


Acceleration 
Spectral 
Density at 
Fundamental 
Mode in g^/ 
Hertz 


Excitation 
Spectral 
Density at 
Fundamental 
Mode in 
(psi) 2 /rad. 
per sec. 


1 

2 


0. 22309 226E-05 
0.22350181E-05 


3072.8148 

3078.4558 


50.011887 

50.103698 


0.14171840E-04 

0.14171840E-04 



Item 

SMS 

Displacement 
in inch 

RMS 

Stress in 
psi 

RMS 

Acceleration 
in g 

RMS 

Excitation 
m psi 

1 

2 

0. 8693253 6E-02 
0.87 852818E-02 

322.63266 

326.04810 

39.575629 

39.582969 

0.090482687 

0.090482687 

b. 

Location x = 

= 11.87 Inches 



y == 

: ^ = 14.594 Inches 



Run No. 4519 Frequency Increment 

n = 33 


Item 

Displacement 
Spectral 
Density at 
Fundamental 
Mode in 
inch 2 /Hertz 

Stress Spectral 
Density at 
Fundamental 
Mode in 
(psi) 2 Hertz 

Acceleration 
Spectral 
Density at 
Fundamental 
Mode in 
g /Hertz 

Excitation 
Spectral 
Density at 
Fundamental 
Mode in 
(psi) 2 /rad 
per sec. 

1 

2 

0 . 56087514E-06 
0.55967865E-06 

938.52039 

936.51828 

125.73464 

125.46641 

0. 14171840E-04 
0.14171840E-04 

Item 

RMS 

Displacement 
in inch 

RMS 
Stress 
in psi 

RMS 

Acceleration 
in g 

RMS 

Excitation 
in psi 

1 

2 

0.0045961768 

0.00453378169 

188.01195 

185.62468 

25.549804 

25.572445 

0.090482687 

0.090482687 


Note: For frequency increment n = 33, the peak of the response spectrum 

occurs at 215*2 Hertz. 


c. 


Location x 

y 

Run No. 6520 


= y ~ 23.75 Inches 

= ^ = 14.595 Inches 

Frequency Increment n = 50 
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Item 


Displacement 
Spectral 
Density at 
Fundamental 
Mode in 


Stress Spectral 
Density at 
Fundamental „ 
Mode in (psi) / 
Hertz 


Acceleration 
Spectral 
Density at 
Fundamental 

n 

Mode in g /Hertz 


Excitation 
Spectral 
Density at 
Fundamental 
Mode in (psi) / 


1 

0.11141607E-05 

1644.1638 

26.010334 

0.14039524E-04 

2 

0. 11142897E— 05 

1644.3542 

26.013345 

0.14039524E-04 


RMS 

RMS 

RMS 

RMS 

Item 

Displacement 

Stress in 

Acceleration 

Excitation 


in inch 

psi 

in g 

in psi 

1 

0.0062653206 

240.68110 

28.105585 

0.090482687 

2 

0.0062120453 

238.63453 

28.061311 

0.090482687 

d. 

Location: x = 

^ = 11.87 Inches 




y = 

£ = 29.187 Inches 

2 




Run No. 7412 

Frequency Increment n = 50 



Di sp lacement 

Stress Spectral 

Acceleration 

Excitation 


Spectral 

Density At 

Spectral 

Spectral 

Item 

Density at 

Fundamental 

Density at 

Density at 


Fundamental 

Mode in 

Fundamental 

Fundamental 


Mode in 

(psi) ^/Hertz 

Mode in 

Mode in 


inch^ /Hertz 


g^/Hertz 

(psi) ^ /rad 





per sec. 

1 

0.11158915E-05 

1670.4033 

26.050740 

0. 14039524E-04 

2 

0 . 11161287E-05 

1670.7584 

26.056278 

0.14039524E-04 


RMS 

RMS 

RMS 

RMS 


Displacement 

Stress 

Acceleration 

Excitation 

Item 

in inch 

in psi 

in g 

in psi 

1 

0.0063428457 

245.40528 

36.238683 

0.090482687 

2 

0.0064173551 

248.28806 

36.224862 

0.090482687 


Note: For frequency increment 

ocurrs at 217.4 Hertz 


n = 50, the peak of the response spectrum 
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3.4 Average Responses Over the Whole Structure 


The formulations for the average responses over the whole structure are 
given in Section 2.9. The 4-series programs, RSRPC4, RFRPC4, and RSFRP4, are 
written to calculate the average responses over the complete structure as well 
as those at any local point. The input data for these programs are the same 
as those for the 1-series program and Program RANDOM (see Section 5.7). 

Following are the results of typical runs of Program RSRPC4 with input data des- 
cribed in Section 5.7. Data in item 1 are the local response at local point of the 
panel and those in item 2 are the average responses over the whole panel. Typical 
plots of average response spectral densities are shown in Figures 35 through 
46, with plots of responses at a local point for comparison. 

% 

a. Local point: x = — = 23.75 Inches 

y = k = 29.187 Inches 

y 2 

Run No. 8009 Frequency Increment n = 33 

Radius of curvature a = 100 Inches 
Fundamental mode Frequency 167.7 Hertz 


Displacement 

Spectral 

Item Density at 

Fundamental 
Mode in 
inc}i 2 /Hertz 


Stress Spectral 
Density at 
Fundamental 
Mode i* 1 
(psi) 2 /Hertz 


Acceleration 
Spectral 
Density at 
Fundamental 
Mode in 
g 2 /Hertz 


Excitation 
Spectral 
Density at 
Fundamental 
Mode in 
(psi) 2 /rad. 
per sec. 


1 


0.93176196E-05 

12833.847 

77.012151 

0.16478916E-04 

2 


0. 23318242E-05 

3211.794 

19.27303 

0.16478916E-04 

Ratio 

of 

0.25 

0.25 

0.25 


Item 2 
1 

to 

RMS 

RMS 

RMS 

RMS 



Displacement 

Stress in 

Acceleration 

Excitation 

Item 


in inch 

psi 

in g 

in psi 

1 


0.15389093E-01 

571.13519 

43.193045 

0.090482687 

2 


0. 78098409E-02 

289.84652 

25.131164 

0.090482687 

Ratio 

of 

0.507 

0.507 

0.58 



Item 2 
to 1 
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For this particular case, it is seen that the root-mean-square responses 
at the center of the panel are approximately twice the average responses over 
the whole panel. 

b. Local Point: x = ■— - 23.75 inches 

y = — = 14.595 inches 

4 

Run No. 6072 Frequency Increment n = 33 

Radius of curvature a = 70 Inches 

Fundamental Mode Frequency : 215.2 Hertz 


Item 

Displacement 
Spectral 
Density at 
Fundamental 
Mode in 
inch 2 /Hertz 

Stress Spectral 
Density at 
Fundamental 
Mode in 
(psi) 2 /Hertz 

Acceleration 
Spectral 
Density at 
Fundamental 
. Mode in 
g 2 /Hertz 

Excitation 
Spectral 
Density at 
Fundamental 
Mode in 
(psi) 2 /rad .per 

1 

0. 11195455E-05 

1652.1102 

25.097502 

0.14171840E-04 

2 

0.55951962E-06 

825.23883 

12.536351 

0. 14171840E-04 

Ratio 

0.501 

0.501 

0.501 



EMS 

EMS 

EMS 

EMS 

Item 

Displacement 

Stress in 

Acceleration 

Excitation 


in inch 

psi 

in g 

in psi 

1 

0.62654106E-02 

240.68455 

28.105819 

0.090482687 

2 

0.44660961E-02 

171.56422 

23.146045 

0.090482687 

Ratio 

0.712 

0.712 

0.822 


c. 

£ 

Local Point: x = 

= 11.87 inches 




li 

= 29.187 inches 




Run No. 6654 Frequency Increment 

n - 33 



Radius of Curvature a = 70 inches 
Fundamental Mode Frequency: 215.2 Hertz 


sec. 
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Item 

Displacement 
Spectral 
Density at 
Fundamental 
Mode in 
inch 2 /Hertz 

Stress Spectral 
Density at 
Fundamental Mode 
in (psi) 2/Hertz 

Acceleration 
Spectral 
Density at 
Fundamental 
Mode in g 2/ 
Hertz 

Excitation 
Spectral 
Density at 
Fundamental 
Mode in (psi) 2 
rad. per sec. 

1 

0.11169505E-05 

1671.9886 

25.039328 

0.14171840E-04 

2 

0.55961962E-06 

837.10857 

12.536351 

0.14171840E-04 

Ratio 

0.501 

0.501 

0.501 



RMS 

RMS 

RMS 

RMS 


Displacement 

Stress 

Acceleration 

Excitation 

Item 

in inch 

in psi 

in g 

in psi 

1 

0. 63429322E-02 

245.40861 

36.238836 

0.090482687 

2 

0. 4460961E-02 

172.79365 

23.146045 

0.090482687 

Ratio 

0.704 

0.704 

0.639 



Comparisons of the average responses over the whole panel and the responses 
at several local points for several panel configurations and subjected to 
different spectra of excitation pressure have been investigated. The results 
can be summarized as follows: 

. The ratio of the average response over the whole panel to the 

local response at a certain point is affected by many factors - 
the panel configuration and properties, the excitation spectrum and 
the position of the local point, etc. 

. Because of the relations in the formulations, the ratios of the 

average response spectral densities to the local ones at a certain 
point for any frequency should be equal. This can serve as a check 
for the computation. 

. The ratio of the average displacement and stress overall root-mean- 

square responses to the local ones at a certain point are approximately 
equal to the square root of the ratio of the average displacement and 
stress spectral densities at fundamental mode to the local ones. This 
is an indication that the fundamental mode dominates the overall 
responses. 

. Because of the assumption of linear structure, the ratios of the 
average displacement and stress to the local ones at a certain 
point are always the same. 
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3.5 Estimation By Formulas of Single Degree of Freedom 


3.5.1 Derivation of Formulas 

Since the fundamental mode contributes the most of the response, 
the responses of a continuous structure subjected to random loading may be 
estimated by assuming the structure as a single-degree-of-freedom system 
vibrating at its fundamental frequency. For a single degree of freedom system, 
it is well known that the following relation holds for the spectral densities 
of the displacement response and the excitation force: 


W«> = % p (a)|H(u ))| 2 ' 


(3.5-1) 


where 


§ 

ww 




displacement spectral density 


*pp<“> 


excitation spectral density 


The square of the magnitude of the frequency response function is given by 


m 2 f(cof-m 2 ) 2 + (2Cm 0 u) 2 ] 

where 

m = mass of the system 
<*> 0 = frequency of the system 

to = independent frequency variable 

? = damping ratio of the system 



By definition, the mean-square displacement is given by 


2 

w => 


/ 4> (to) do> 

J ww 


= / * (w)|h(w)| dco 

o " 


(3.5-2) 


(3.5-3) 
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When the excitation spectral density is constant for all frequency, the inte- 
gration in Equation (3.5-3) can be carried to yield the mean-square displace- 
ment as 


w 


2 


ir§. 


£EL 




(3.5-4) 


When the excitation spectral density is not constant, Equation (5.3-4) still 
gives very good approximation of the mean-square displacements. 


The half-power bandwidth at the fundamental mode is 


B c = 2 4<o 0 


(3.5-5) 


The displacement spectral density at resonance can be obtained by setting 
m = m 0 in Equation (3.5-2) and substituting the result into Equation (3.5-1) 
as 


*ww (w o> 




(2£u)2mV 


(3.5-6) 


The mean-square displacement may also be approximated by the product of the 
displacement spectral density and the bandwidth at resonant frequency: 


w 


= B 0 $ w (m 0 ) 


(3.5-7) 


The acceleration spectral density of the single degree of freedom system will 
be given by 


$ ww(“) “ (<*0 = “ 4$ pp ( w ) |h(“) i 2 


(3.5-8) 


where 


^ww( w ) is given by Equation (3.5-1) . 

The acceleration spectral density at resonance can be obtained by setting 
w = m 0 in Equation (3.5-8); 


%wK) 


a 4 

o 


w - ^ww (%) 


~PP 

(2?m) 2 


(3.5-9) 
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The mean-square acceleration is given by 


G 2 = / 


(3.5-10) 


In general, the excitation spectral density ®pp(oi) can not be .expressed in an 
analytical form. Hence, the integration in Equation (3.5-10) has to be evalu- 
ated numerically. However, approximation to the mean-square acceleration can 
be obtained by miltiplying the spectral density at resonance with the half- 
power bandwidth at resonance frequency: 


B $—(o)) 
o ww v ' 


tin 


*pp(<V) 

2?m 2 


4 2 

= “o w 


(3.5-11) 


When these equations are used to estimate the responses of a continuous 
structure, the fundamental frequency will be used for m 0 , the spectral density 
of the excitation pressure will be used for $pp, the mass per unit area of the 
structure will be used for m, and the damping ratio for the fundamental mode 
will be used for £. 

If only a portion of the structure is subjected to the excitation, the 
square of the ratio of the area under excitation to the total area may be used 
to multiply the results to give the responses: 


$ 

ww 


(to) = 


w 


2 


w 


2 


^ww^o^ = 


( f-)%>p (<0) l H(u) 


(«■; 


^PP^) |H(«) | 2 dto 


(f) 


r \ 2 ir$. 


PP 


3 ? 


(r) 


•\2 * pp ( “ 0 > 


(2Cm o 2 m) 2 


(3.5-12) 


(3.5-13) 


(3.5-14) 


(3.5-15) 
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(3.5-16) 


$ (w) 

ww'' J 



*pp<“o> 


<2SJn)‘ 



u o. ^ww 


(“n) 


G 


2 



“o i ppK>> ,/ Slf m 4»2 
2 C m 2 . K S 1 ° 


(3.5-17) 


3.5.2 Example 

As an example to illustrate the application, the above formulas 
are used to estimate the responses of the panel described in Section 3.4. In 
this case: 

The spectral density of the excitation pressure corresponding to the 
fundamental frequency is 


,-5 


PP 


= 1.648 x 10 (psi) /rad. P er sec. 


The fundamental frequency is 

io 0 = 167.7 Hertz = 1055 rad/sec. 

The damping ratio is 

e = 0.04 


The smeared-out mass per square inch is 


m 


ph + p* h’ 


2.51 x 10 4 x 0.05 + 2.60 x 10' 4 x 0.05 
2.555 x 10“ ^ lbf-sec 2 /in^ 


The area of the panel is 


S = = 58.375 x 47.5 = 2871 inch 2 


The area of the panel exposed to the excitation pressure is 
S T = b f £’ = 50 x 45 = 2250 inch 2 
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Substitution of the above values into Equation (3.5-14) gives the mean-square 
displacement as 


2 _ tt( 1.648 x 10 5 ) . / 2250 \ 2 

4(0.04) (1054) 3 (2.555 x 10~ 5 ) Z \ 2871 / 

= 2,6 x 10“ 4 inch 2 

or the root-mean-square displacement is 

w = (2.6 x 10 4 ) = 0.0161 inch 


This agrees very closely with the actual root-mean-square displacement at the 
center of the panel calculated by Program SSEPC4 which is 


w ** 0.01539 inch 


Substitution of the above values into Equation (3.5-6) gives the displacement 
spectral density at the fundamental frequency as 


$ 

ww 


C«o) = 


1.648 x 10 

(2 x 0.04 x 1055 2 x 2.555 x 10T 5 ) 2 



2 


= 1.96 x 10“ 6 


in . / rad. 


per sec. 


12.3 in?/Hertz 


The actual displacement spectral density at the fundamental frequency calculated 
by Program RSKPC4 is 

$ ww^“cj) = 9 *3176 in. 2 /Hertz 


From Equation (3.5-5) the bandwidth of the fundamental frequency is 


B q = 2 x 0.04 x 1055 = 84.4 rad. /sec. 
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Substituting these values into Equation (3.5-7) gives mean-square displacement as 


w 2 = 84.4 x 1.96 x 10 6 

= 1.655 x 10 4 inch^ 

Hence, the root-mean-square displacement is 
w = (1.655 x 10 _4 ) 1/2 

= 0.01286 inch 

From Equation (3.5-9), the acceleration spectral density at the fundamental 
frequency is 

*55 (%) = 1054 4 (1.96 x 10" 6 ) = 16.3 g 2 /rad. per sec. 


From Equation (3 .5-11) j the mean-square acceleration is 

G 2 = 84.4 X 16.3 = 1377 g 2 

Hence, the root-mean-square acceleration is 
G = 37.05 g 

The actual root-mean-square acceleration calculated with Program RSRPC4 is 
43.193 g. 
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3.6 Number of Modes and Modal Density 


Analytical expressions for the total number of modes and the modal density 
are involved. Only approximate solutions to simple cases can be found in the 
literature. When the natural frequency of a structure have been calculated from 
the frequency equations, the total number of modes and modal density can be 
solved easily with the computer. 

Let f be the frequency of a two dimensional structure as solved from the 
frequency equations or some other method. Arrange fj^ in the ascending order in 
a single array. 


f 


n 


n 1, 2 , ... 


The curve of n versus f will give the total number of modes up to any fre- 
quency range. 

The slope of the n-f curve will yield the modal density as a function of 
frequency. In the calculation, the slope of the curve is approximated by 


dn _ 
df 

f = f. 
J 



where 


(3.6-1) 


4nj = a, + i - nj - 1 

i£ J - fj+i-fj <3-6-2> 

j = 1,2,3... 

The plot of Anj/Afj versus fj gives the variation of modal density as a function 
of frequency. 


Program NFU0P1 is written to investigate the total number of inodes and the 
modal density. Figures 41 and 42 are typical plots of the total number of modes 
and the modal density. These curves are very useful in the analysis of the 
dynamic characteristics of the structure. It is seen that the plot of the 
total number of mode is a rather smooth curve while the modal density curve is 
not. The peaks in the modal density plot indicate where two frequencies are very 
close together. Sudden change in plot of the number of modes will imply some 
peculiar dynamic characteristics of the structure. 
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3 . 7 Joint Acceptance and Correlation Coefficients 


Because of their influence on the responses and the complexity of the 
analytical expressions for the joint acceptance and the cross spectral density 
of the generalized force, Program JAESR1 is specially written to investigate 
the behavior of these two quantities as a function of the frequency with mode 
indices as parameter for both correlation functions described in Section 2.3. 

This program also computes and plots the curves of the correlation coefficient 
as a function of separation distance with the frequency as parameter-. Figures 
43 and 44 are typical curves of the exponentially decaying correlation coefficient 
(Equation 2.3-3). Representative curves of the sinusoidal decaying correlation 
coefficient (Equation 2.3-4) are shown in Figures 45 and 46. Figures 47 through 
52 are typical curves of these quantities for the exponentially decaying corre- 
lation. Typical plots of the two quantities for the sinusoidal decaying corre- 
lation are shown in Figures 53 through 58. It is noted that is identical 

with IjlcjR • The curves of these two quantities are more regular for the ex- 
ponentially decaying correlation than for the sinusoidal decaying correlation. 

3.8 Vibro-Acoustic Transfer Function 


The vibro-acoustic transfer function at any point of a structure is defined 
as the ratio of the response to the excitation. When the response and the ex- 
citation are expressed in decibel scale, the vibro-acoustic transfer function in 
decibel scale at any point r will be the excitation in decibels minus the response 
in decibels : 


where 


T”f,f) 


s " g (r,f) 



(3.8-1) 





£,f> 


Vibro-acoustic transfer function in decibels with 
excitation pressure referenced 0,00002 Newton/meter^ 
and acceleration response referenced g 

Acceleration spectral density in decibels referenced 
g given by Equation (2.7-2b) in Section 2.7 


Spp(f) = The spatial homogeneous excitation spectral density 
1 ' in (psi) ^/Hertz given by Equation (2.6-3) in Section 

2 . 6 . 


Program RSRPC2 is written to investigate the vibro-acoustic transfer' function 
of four edges simply-supported rectangular cylindrical shell panels cross re- 
inforced with ribs and stringers. Typical plots of this program are shown in 
Figures 59 and 60. 
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IV. COMPARISON OF COMPUTED RESULTS WITH TEST DATA 


The developed computer programs have been used in the research project 
which Chrysler Huntsville Operations was currently conducting for Marshall Space 
Flight Center under Contract No. NAS 8-21425. The object of this project is to 
develop comparative analysis of acoustic testing techniques and to determine 
the time for acoustic qualification test at level other than specified levels. 
Comparison of the computed results with test data shows good agreement. 

One of the test specimens in Project NAS8-21425 was a 4 in. x 13 in. } 

0.2 in. thick aluminum flat plate. The plate is knife-edges supported on 
four sides and is subjected to high level acoustic pressure until fatigue fail- 
ure occurs . The input spectrum of the acoustic pressure is shown in Figure 
61 . The following input data are used in Program RFRPC1 to calculate the 
responses at six locations of the plate. Figure 62 shows the calculated 
acceleration spectral density compared with test data at location number 9 , the 
coordinates of which are x = 6.5 inches and y - 1 inch. 


Notation 

In 

Formulas 

E 

E' 


’jk 

a 

a 

b 

V 


Definition 

Fortran and 

Notation Description 

E Young's modulus of panel skin 

EP Young's modulus of stiffeners 

VIP Poisson's ratio of panel skin 

RHO Mass density of panel skin 

RHOP Mass density of stiffeners 

Cl Damping ratio of panel 

RAD Radius of panel 

PL Axial length of panel 

B Width of panel 

PLP Length of panel subjected to 

excitation 


Data Used 
In 

Test Run 

10 x 10 6 lbf/in. 2 
0 

0.3 

2.58 x 10 - ^ lbf-sec 2 /in.^ 
0 

0.065 
inches 
13 inches 
4 inches 
13 inches 


b' 


a l 


BP 

Width of panel subjected to 
excitation 

4 inches 

AL1 

Spacing of width-direction 
stiffeners 

®> inches 


60 



Notation 

In 

Formulas 


Fortran 

Notation 


Definition 

and 

Description 


BL1 Spacing of length-direction 
stiffeners 


h 

h' 


HS Thickness of panel skin 

HP Smeared- out thickness of 

stiffeners 


H2 Largest height of stiffeners 

at point investigated 

All Moment of inertia of one length- 
direction stiffener with respect 
to neutral axis 

AI2 Moment of inertia of one width- 
direction stiffener with respect 
to neutral axis 


y 

a. 


X Coordinate of r 

Y Coordinate of r 

A1 Correlation decay constant 

in axial length-direction 

A2 Correlation decay constant 

in circumferential width- 
direction 


s 3 t (f) 


C Speed of sound in air 

S3RD One-third octave pressure 

level spectrum of excitation 


n 


FINN 


One-nth octave increment 


Data Used 
In 

Test Run 

03 inches 

0 . 2 inch 
0 inch 

0 

0 inch^ 

0 inch^ 

6 . 5 inches 

1 inch 
10 

0 

13500 in. /sec. 
See Figure 61 

33 
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V. DESCRIPTION OF THE COMPUTER PROGRAMS 


5.1 Program RANDOM 

Program RANDOM is a combination of the three main programs RSRPCl, RFRPC1, 
and RSFRP1 described in the following sections. With this .program a single loading 
of the input data will be sufficient to obtain the responses of a rectangular 
cylindrical shell panel cross-reinforced with stiffeners and subjected to three 
different boundary conditions : all edges simply- supported , all edges clamped, 

and two opposite edges simply-supported with other two clamped. Controls are 
also provided to run any one of the three boundary conditions . 

The input data for this program are (a) the material constants, the geometric 
dimensions and properties of the shell panel and stiffeners, (b) the one- third 
octave spectrum of the excitation pressure, (c) the x and y coordinates of the 
point of interest and, (d) some control constants. (see Section 5.7). 

The output data of this program include: 

a. All the input data with nomenclature. 

b. The input one-third octave spectrum of the excitation pressure, this 
excitation spectrum is also converted into spectral density both in decibel/ 

Hertz and (psi)2/Hertz, tabulated and plotted in graphical forms. 

c. The natural frequencies of the panel both in Hertz and radians/second. 

d. The spectral densities of displacement in (inch) /Hertz, of stress 
response in (psi)^/Hertz, of acceleration response in g^/Hertz and the excitation 
spectral density in (psi)2/rad. per sec. In addition to tabulation, all three 
response spectral densities are plotted out in graphical forms. 

e. The mean-square and the root-mean-square values of the responses 
and excitation. 

9 

f. The constants Q x , Qy, and Qw (QX, QY and QW) and the constant y 
(GAMMA2) to convert the displacement spectral density into stress spectral 
density. 


g. Some values of the joint acceptance squared of all combination of 
modes for the beginning frequency. 

5.2 Programs RSRPCl, RSRPC2, RSSPC3 and RSRPC4 

These four programs are writen and designed to calculate the responses 
of the simply-supported cylindrical shell rectangular panel cross-reinforced 
with ribs and stringers for various specific purposes. The input data for these 
programs are the same as for Program RANDOM described in Section 5.7. The 
output data for these four programs are the same as Program RANDOM but with 
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extra output for each individual program. 

Program RSRPC1 calculates the spectral densities of the displacement 
response in (inch) 2/Hertz, the stress response in (psi) 2/Hertz , and the accel- 
eration response in g2/Hertz. Therefore the output data are identical with 
Program RANDOM run in the case of all edges simply-supported. 

In addition to those as in Program RSSPC1, the output data of Program 
RSRPC2 include the acceleration spectral density in decibels referenced gravity 
acceleration (g) and the vibro-acoustic transfer function as a function of fre- 
quency. Both the transfer function and the decibel scale acceleration spectral 
density are also plotted out in graphical forms. 

Program RSRPC3 is designed to investigate the contribution of the cross 
terms to the vibrational responses. Both the responses of all terms summation 
and cross terms neglected are tabulated and plotted for comparison. Formulations 
for responses with all terms summation and cross terms neglected and analysis 
of the results are given in Section 3.3. Typical plots of this program are' shown 
in Figures U through 34. 

Program RSRPC4 is a modification of Program RSRPC1 and is written to calcu- 
late responses at any point of the structure and the average responses over the 
whole structure. Formulations for the average displacement response over the panel 
are derived in Section 2.9. Analysis of the results is given in Section 3.4. 
Typical plots of this program are shown in Figures 35 through 40. 

5.3 Program RFRPC1 and RFRPC4 

These programs are written to calculate the vibrational responses of the 
four edges clamped cylindrical shell rectangular panel cross-reinforced with 
ribs and stringers. The input data for these programs are the same as those for 
Program RANDOM described in Section 5.7. The output data for these programs are 
essentially the same as Program RANDOM except Program RFRPC4 also gives the 
average responses over the whole panel. 

Program RFRPC1 calculates and plots the spectral densities of the displace- 
ment response in (inch) 2/Hertz, the stress response in (psi)2/Hertz, and the 
acceleration response in g^/Hertz. Therefore, the output data are identical 
with Program RANDOM run in the case of all edge clamped, except the frequency 
range can be higher and the frequency increment can be smaller. 

Program RFRPC4 is a modification of Program RFRPC1. In addition to the 
output of Program RFRPC1, Program RFRPC4 calculates the average responses over 
the whole structure. Formulations for the average responses over the whole 
panel are derived in Section 2.9. 
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5 . 4 Programs RSFRPl and RSF5P4 


These programs are written to calculate the vibrational response of the 
two opposite edges simply-supported while other two clamped rectangular cylin- 
drical shell panel cross-reinforced with ribs and stringers. The input data 
for these programs are the same as Program RANDOM described in Section 5.7. The 
output data for these programs are the same as Program RANDOM except that Program 
RSFHP4 also calculates the average responses over the whole panel. 

Program RSFRPl calculates and plots the spectral densities of the displace- 
ment response in (inch) 2/Hertz , the stress response in (psi) 2/Hertz, and the 
acceleration response in g^/Hertz. Therefore, the output data are identical 
with Program RANDOM run in the case of two opposite edges simply-supported while 
other two clamped, except with higher frequency range and smaller frequency 
increment . 

Program RSFRP4 is a modification of RSFRPl. In addition to calculating the 
local responses at any point of the structure, Program RSFRP4 also calculates and 
plots the average responses over the whole panel. Formulations for the average 
responses over the whole panel are derived in Section 2.9. 


5.5 Program JARSR1 

This is a computer program written to study the correlation coefficient as 
a function of separation distance, the joint acceptance and the normalized cross 
spectral density of the generalized force as a function of frequency for various 
correlation, functions. The correlation functions investigated and the derivation 
of the expressions for the joint acceptance and *the normalized cross spectral 
density of the generalized force are given in Sections 2.3 and 2.4. Typical 
plots of the correlation coefficients, the joint acceptance and the normalized 
'cross spectral density of the generalized force are shown in Figures 43 through 
58. 


5.6 Program NFU0P1 

Program NFU0P1 is to calculate the total number of modes and the modal 
density. The input data to this program are the material constants, the geo- 
metric dimensions and properties of the structure required in the frequency 
equations. The output data of this program are the tabulated and plotted number 
of modes and the modal density as a function of frequency. Formulations of this 
program are given in Section 3.6. Sample results are given in Figures 41 and 
42. The analysis of the number of modes and modal density is very useful 
in the evaluation of the frequency equations. It enables the structures engineer 
to determine the accuracy and the behavior of the frequency equations. This 
program can be modified to apply to any structure. 
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5.7 The Input Data 


The input data for these programs include: (a) the material constants, the 

geometric dimensions and properties of the shell panel and the stiffeners, (b) 
the one-third octave pressure level spectrum and the spatial correlation properties 
of the excitation pressure field, (c) the coordinates of the point concerned, and 
(d) some control constants. These are listed as follows: 


Notation 


Definition 

Typical Data Used 

In 

Fortran 

And 

In 

Formulas 

Notation 

Descritption 

Test Run 

E 

E 

Young’s modulus of, panel skin 

10 x 10 6 lbf/.in. 2 

E’ 

EP 

Young's modulus of stiffeners 

12 x 10 6 lbf/in. 2 

V 

VIP 

Poisson’s ratio of panel skin 

0.3 

P 

EHO 

Mass density of panel skin 

2.51 x 10 4 lbf-sec' 

P f 

KHOP 

Mass density of stiffeners 

‘ 2.60 x 10~ 4 lbf-sec' 


Cl 

Damping ratio of panel’ 

0.04 

a 

RAD 

Radius of panel 

100 inches 

£ 

PL 

Axial length of panel (along 
x-axis) 

47.5 inches 

b 

B 

Width of panel (along y-axis) 

58.375 inches 

£' 

PLP 

Length of panel subjected to 
excitation 

45 inches 

b’ 

BP 

Width of panel subjected to 

50 inches 


- 

excitation 


a l 

AL1 

Spacing of width-direction . 
stiffeners 

11.9 inches 

b l 

BL1 

Spacing of length-direction 
stiffeners 

14.6 inches 

h 

HS 

Thickness of panel skin 

0.1 inch 

h' 

HP 

Smeared-out thickness of 
stiffeners 

0.08 inch 

h 2 

H2 

Largest height of stiffeners 
at point investigated 

1.0 inch 
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Notation 

In Fortran 

Formulas Notation 


Definition 

and 

Description 


Data Used 
In 

Test Run 


I 


1 


I 


2 


x 


All Moment of inertia of one length- 

direction stiffener with respect 
to neutral axis 

AI2 Moment of inertia of one width- 

direction stiffener with respect 
to neutral axis . 

X Coordinate of r 


0.5 


inch 


4 


0 . 6 inch^ 


23 . 75 inches 


y 


-y 

Y Coordinate of r 


28.1875 inches 


% 


A1 Correlation decay constant 

in axial length-direction 


10 


A 2 


c 


A2 Correlation decay constant 

in circumferential width- 
direction 

C Speed of sound in air 


0 


13500 in. /sec 


s 3r ( f ) 


S3RD One-third octave pressure 

level spectrum of excitation 


See Figures 63 
and 64 and the 
following table 


f F3RD Input spectrum frequencies 

n FINN One-nth octave increment 


33 


Note: Refer to Figure 1 for identification of geometrical dimensions. 
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The one-thrid octave level spectrum of the excitation pressure input data 
in a typical test run is as follows: 


0ne-3RD Octave 
Mean Frequency 


F3RD 


0ne-3RD Octave Level 
In Decibels 
S3rd 


5.000 
6.300 

8.000 
10.000 

12.500 
16.000 
20.000 

25.000 

31.500 

40.000 

50.000 

63.000 

80.000 

100.00 

125.000 

160.000 

200.000 

250.000 

315.000 

400.000 

500.000 

630.000 

800.000 

1000.000 

1250.000 

1600.000 
2000.000 

2500.000 

3150.000 

4000.000 

5000.000 


129.000 

131.000 

133.000 

135.000 

136.000 

138.000 

139.000 

140.500 

142.000 

143.000 

144.000 

145.000 

145.500 

145.000 

146.500 

146.500 

147.000 

147.000 

146.000 

145.500 

145.500 

144.000 

142.500 

141.000 

139.500 

137.500 

136.000 

134.500 

133.000 

131.000 

129.500 


Note: Point frequencies can be any discrete frequencies with any increment. 

The above spectrum is the preliminary acoustic test specification for 
components of Saturn IB vehicle in Subzone 7-1 for static firing, (See reference 
6, Page 64). 
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The input data spectrum of the excitation pressure is converted into 
spectral density in decibels/Hertz , (psi)^/Hertz and (psi) ^/rad . /sec, by the 
computer and plotted as in Figures 63 and 64 . 


Frequency 

In 

Hertz 


Spectral Density Spectral Density 

In Decibels /Hertz In PSI SQ/RAD/SEC 

SPPF FIPW 


5 

6.3 

8 

10 

12.5 

16 

20 

25 

31 

40 

50 

63 

80 

100 

125 

160 

200 

250 

315 

400 

500 

630 

800 

1000 

1250 

1600 

2000 

2500 

3150 

4000 

5000 


128.365 

0.95653E-05 

129.361 

0. 12032E-04 

130.324 

0.15017E-04 

131.354 

0. 19040E-04 

131.885 

0.21516E-04 

132.313 

0.23744E-04 

132.344 

0.23913E-04 

132.875 

0.27023E-04 

133.371 

0.30294E-04 

133.334 

0.30034E-04 

133.365 

0.30248E-04 

133.361 

0.30222E-04 

132.824 

0.26704E-04 

132.354 

0.23970E-04 

131.885 

0.21516E-04 

130.813 

0.16809E-04 

130.344 

0.15088E-04 

129.375 

0.12071E-04 

127.371 

0 . 76095E-05 

125.834 

0.53408E-05 

124.865 

0.42727E-05 

122.361 

0.24006E-05 

119.824 

0.13384E-05 

117.354 

0.75800E-06 

114.885 

0.42930E-06 

111.813 

0.21162E-06 

109.344 

0.119 85 E-06 

106.875 

0. 67878E-07 

104.371 

0.38138E-07 

101.334 

0.18950E-07 

98.865 

0. 10732E-07 
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5.8 The Output Data 

The output data of the programs for the responses are as follows: 

a. All the input data. 

b . The excitation spectrum expressed in one- third octave level in 
decibels, in spectral density in decibels per Hertz, and in (psi)2 per radian 
per second. 


c. The natural frequencies of the structure both in Hertz and radians 
per second. 


d. The displacement, the stress, the acceleration and the excitation 
pressure spectral densities both in tabulations and plots. 

e. The mean-square and the root-mean-square values of the responses 
and the excitation. 

f. Some values of the joint acceptance. 

g. Other special output and plots for each individual program.' 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


6.1 Conclusions 

6.1.1 In the calculation of the random structural vibrational responses 
due to the fluctuating pressure environments by the normal mode approach, the 
important factors that have to be determined are: 

a. The normal modes 

b. The natural frequencies 

c. The joint acceptance or the cross spectral density of 
the generalized force which in turn depends on: 

(1) The normal modes 

(2) The correlation properties of the pressure field. 

(3) The spectral density of the excitation pressure. 

6.1.2 For simple structures, the determination of the mode shapes and 
frequencies can be obtained to a certain degree of accuracy. For complex 
structures, the mode shapes and frequencies are difficult to determine in 
general. 

6.1.3 Except for simple pressure fields, the correlation properties 
are not easy to determine. Some experiments have been done in this field, 
but scattering of the test data is tremendous. 

6.1.4 Any discrepancy of the factors mentioned in 6.1.1 will affect the 
accuracy of the computed results . 

6.1.5 To develop computer programs for the calculation of this kind of 
random vibrations, the important tasks are to search or develop the necessary 
frequency equations, the normal modes, the analytical expressions for the 
correlation properties and the joint acceptance. The degree of accuracy of 
the expressions for these quantities must be consistent. When these quantities 
are incorporated in the formulations of the programs, the requirement of com- 
puter time should be moderate as to make them practical for use. Such tasks 
have been done in this project. 

6.1.6 One of the important features of the computer programs developed 
here is that any shape of the spectrum of the excitation pressure can be input 
into the programs. Thus, the simulation of the excitation will be as accurate 
as the spectral analysis of the random pressure. 

6.1.7 The discrepancy in the determination of the natural frequencies 
will affect the response spectrum to some degree while it will not affect the 
overall mean-square response very much. 
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6.1.8 It is well known that the determination of the damping properties 
of structures is very difficult. By utilization of the developed computer 
programs and the results of test data, the damping properties of structures 
can be determined. This is another application of these programs. 

6.1.9 In the design stage of structures of aerospace vehicles, desk 
calculation to estimate the responses is very important. The developed 
computer programs are good guides for the designer as to the degree of approxi- 
mation that can be obtained by assuming the structure vibrating at its funda- 
mental mode and estimating the responses by intelligent use of simple formulas. 

6 . 2 Recommendations for Further Investigation 

A study of the random vibration of structures is a very complicated pro- 
ject. There are many additional areas not studied in detail during the per- 
formance of this project that represent opportunities to extend the field. It 
is believed the following recommendations are worthy of further investigations. 

6.2.1 When a shell structure is reinforced by ribs and stringers, in 
addition to the vibrations of the complete structure as a system, there are 
local vibrations of the portions between the stiffeners. Thus the resultant 
vibration of any point should be the superposition of the system vibrations 
and the local vibrations. The computer programs developed in this project 

can be easily modified to account for the local vibrations. When this is done, 
the prediction of the structural vibrations will be significantly improved. 

6.2.2 In the development of the computer programs, it is found that the 
analytical expressions for the joint acceptance are the most complicated to 
incorporate into the programs. Further study can be made to de rive more 
analytical expressions of the joint acceptance for some of the available 
correlation functions and structure mode shapes. For the cases analytical 
expressions are not practical, numerical integration should be incorporated. 

6.2.3 For complex structures, only finite element methods will be 
effective for the determination of the mode shapes and f recjuencies . Tremendous 
progress has been made in the stress analysis of structures with the finite 
element method. However, little has been done in structural random vibration 
with this method. It is believed that application of the finite element method 
to random vibrations of structures provides a promising field of research. The 
program developed in this project can be adapted to the finite element method 
without difficulty. 

6.2.4 The vibro-acoustic transfer function is as important to random 
vibrations as the impedance to harmonic vibrations. Experimental determination 
of the vibro-acoustic transfer functions is expensive. The development of the 
2-series programs to calculate the transfer functions at any point of the 
structure subjected to any excitation pressure is promising. Further development 
in this respect is recommended. 
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6.2.5 In practical applications , it may happen that the average responses 
over the whole structure are of interest. Formulations for the average responses 
have been derived for different boundary conditions. Further study in the pro- 
grams is necessary. 
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VII. NOMENCLATURE 


Notation 


4 2 ' 


^jkmn 


B i 


B jkmn 

^jkmn 

Tj _2l/2n_«-l/2n 
n ^ 


D n = 2 1/n -l 


D 

x 

D 

y 

E 

E' 


F jk (r) 

G(J) 

G 2 (r) 

G 2 

1 

H 


H* 

mn 


I 

I 


1 

2 


Description 

Decay constant along x-axis 
Decay constant along y-axis 
See equation (2.4-29) 

Half -power bandwidth of ith mode 
See equation (2.4-29) 

See equation (2.4-29) 

One-nth octave bandwidth constant referred to 
geometric mean frequency 

One-nth octave bandwidth constant referred to 
lower limit 

Rigidity, equation (2.1-4) 

Rigidity, equation (2.1-4) 

Young's modulus of panel skin 
Young's modulus of stiffeners 
Normal mode 

Root-mean-square acceleration in g 
Mean-square acceleration in g^ 

Mean-square acceleration contributed by the ith mode 
Rigidity, equation (2.1-4) 

Frequency response function 
Conjugate of H mn 

Moment of inertia of one length-direction stiffener 
with respect to neutral axis 

Moment of inertia of one width-direction stiffener 
with respect to neutral axis 
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VII. NOMENCLATURE (Continued) 


Notation 


II, I 2 


I , I 
x’ y 


X x> \ 


X * „ 


I jm 


L kn 


i: 

kn 


T xij’ I yij 


I xijk 


L yijk 


T 

jkmn 


^jkmn 


j2. 

jkmn 


K 


M 


M 


jk 


Description 

Integrals, equation (2.4-6.) 

Integrals, equation (2.4-5); also different 
use in equation (2.9-1) 

See equations (2.4-7) and (2.4-7a) 

See equation (2.4-22) 

See equation (2.4-23) 

See equation (2.4-24) 

See equation (2.4-25) 

See equation (2.4-8) 

See equation (2.4-9) 

See equation (2.4-10) 

Cross spectral density of generalized force 

Normalized Cross Spectral Density of generalized 
Force 

Joint acceptance squared 
Wave number 

Excitation overall pressure level in decibels 

Smeared-out mass per unit area 

Modal mass 

See equation (2.8-3) 

See equation (2.8-4) 
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VII. NOMENCLATURE (Continued) 


Notation 

Q w 


s 

s' 


s pp (f) 


s* <f) 
pp 


S wg (r ’ f) 


S ww< r ’ f > 


ww 


(r,f) 


S G a (r ’ f) 

S 3rd 


T w(r.f> 

X j (x) 

*kCy) 


a 


a l 

b 


b* 



c 

e 

f 


Description 
See equation (2.8-5) 

Area of panel 

Area of panel subjected to excitation 
Excitation spectral density in decibels/Hertz 

2 

Excitation spectral density in (psi) /Hertz 

Acceleration spectral density in decibels 
referenced g 

2 

Displacement spectral density in inch /Hertz 

2 

Acceleration spectral density in g /Hertz 

2 

Stress spectral density in (psi) /Hertz 

One-third octave excitation pressure level 
in decibels 

Vibro-acoustic transfer function in decibels 
See equation (2.2-1) 

See equation (2.2-1) 

Radius of shell 

Spacing of width-direction stiffeners 
Circumferential width of panel 
Width of panel subjected to excitation 
Spacing of length-direction stiffeners 
Speed of sound in medium 
Base of natural logarithm 
Frequency in Hertz 
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VII. NOMENCLATURE (Continued) 


Notation 

g 

h 

h* 

= h + I12 
h 2 

i 

% 

V 

n 

p(r) 

Pr 

2 

Pa 

r 

s, s 1 
w(r) 
w 2 (r) 
x,y 

x j’ y k 

V (r l’ r 2 ,ul) 

$pp( w ) 


Description 

Gravity acceleration 

Thickness of panel skin 

Smeared— but thickness of stiffeners 

Height 

-)■ 

Largest height of stiffeners at r 
(see figure 1) 

/-l ; also index 

Mode indices 

Axial length of panel 

Length of panel subjected to excitation 
One-nth octave increment 
Pressure at r 
Referenced pressure 

2 

Overall mean-square pressure in (psi) 

Position vector 
Number of divisions 
Root-mean-square displacement 
Mean-square displacement 
Cartesian coordinates of r 
See equation (2.4-15) 

Cross spectral density 

2 

Excitation spectral density in (psi) /rad. per sec. 
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VII . NOMENCLATURE (Concluded) 


Notation 

W«) 

-V 

%rw^ r > w ) 

*og£> w > 
Oj_ = A-jK£ 

a 2 = A 2 K ^ 
o ■* 

T 2 (r) 



n 

*1 = Ke 

1 2 = Kb 

v 

§ 


p( r l> r 2 ) 

o(?) 


9 

o (r) 


Description 

2 

Displacement spectral density in inch /rad. per 
sec. 

2 

Average displacement spectral density in inch / 
rad. per sec. 

2 2 

Acceleration spectral density in (inch/sec ) / 
rad. per sec. 

2 

Stress spectral density in (psi) /rad. per sec. 
See equation (2.4-13) 

See equation (2.4-13) 

Constant to convert displacement spectral' density 
into stress 

damping ratio 

Separation distance along y-axis 
See equation (2.4-13) 

See equation (2.4-13) 

Poisson’s ratio 

Separation distance along x-axis 
Mass density of panel skin 
Mass density of stiffeners 
Correlation coefficient 

Root-mean-square stress in psi 

mean-square stress in (psi) 2 


(0 

10 jk 


frequency in rad/sec. 

Natural frequency in rad/sec. 


77 



VIII. REFERENCES 


1. Powell, Alan, "On the Fatigue Failure of Structure Due to Vibrations 
Excited by Random Pressure Fields", Journal of Acoustical Society of 
America, Volume 30, No. 12,, December 1958, Pages 1130-1135. 

2. Lin, Y. K. , "Probabilistic Theory of Structural Dynamics" , McGraw- 
Hill, 1967. 

3. Lee, T. N., "Computer Program for Prediction of Structural Vibrations 
to Fluctuating Pressure Environments," Monthly Progress Report 'No. 4, 
Contract No. NAS8-21403, Chrysler Corporation Space Division, Huntsville 
Operations, October 1968. 

4. Lee, T. N., "Computer Program For Prediction of Structural 
Vibrations to Fluctuating Pressure Environments," Monthly Progress 
Report No. 3, Contract No. NAS8-21403. Chrysler Corporation Space 
Division, Huntsville Operations, September 1968'. 

5. Lee, T. N., "Computer Program For Prediction of Structural Vibrations 
To Fluctuating Pressure Environments," Monthly Progress Report No. 11, 
Contract No. NAS8-21403, Chrysler Corporation Space Division, Huntsville 
Operations, May 1969. 

6. Sailors, C. W. and Lehner, F.A. , "Preliminary Acoustic Test Specifications 
For Components on Saturn IB Vehicle", Technical Report HSM-N02-68, 

Chrysler Space Division, 2 June 1968. 


78 



IX. APPENDIX 


9.1 Figures 
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Note: Numbers in parenthesis indicate mode indices and frequencies 

FIGURE 2 ACCELERATION SPECTRAL DENSITY AT CENTER OF FOUR EDGES SIMPLY- 
SUPPORTED CURVED RECTANGULAR PANEL CROSS-REINFORCED WITH 
STIFFENERS 
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Note: Numbers in parenthesis indicate mode indices and frequencies. 

FIGURE 3 ACCELERATION SPECTRAL DENSITY AT CENTER OF FOUR-EDGES CLAMPED 
CURVED RECTANGULAR PANEL CROSS-REINFORCED WITH STIFFENER 
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FIUGRE 4' ACCELERATION SPECTRAL DENSITY AT- CENTER OF TWO OPPOSITE 
EDGES SIMPLY-SUPPORTED AND OTHER TWO CLAMPED CURVED 
RECTANGULAR PANEL CROSS-REINFORCED WITH STIFFENERS 
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FIGURE 5 DISPLACEMENT SPECTRAL DENSITY AT CENTER OF FOUR 
EDGES SIMPLY-SUPPORTED CURVED 'RECTANGULAR PANEL 
CROSS-REINFORCED WITH STIFFENERS 
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FIGURE 6 DISPLACEMENT SPECTRAL DENSITY AT CENTER OF FOUR EDGES 
CLAMPED CURVED RECTANGULAR PANEL CROSS-REINFORCED 
WITH STIFFENERS 
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FIGURE 7 DISPLACEMENT SPECTRAL DENSITY AT CENTER OF TWO OPPOSITE 
EDGES SIMPLY-SUPPORTED AND OTHER TWO CLAMPED CURVED 
RECTANGULAR PANEL CROSS- REINFORCED WITH STIFFENERS 
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FIGURE 8. STRESS SPECTRAL DENSITY AT CENTER OF FOUR EDGES SIMPLY- 
SUPPORTED CURVED RECTANGULAR PANEL CROSS-REINFORCED 
WITH STIFFENERS 
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FIGURE 9. STRESS SPECTRAL DENSITY AT CENTER OF FOUR EDGES CLAMPED 

CURVED RECTANGULAR PANEL CROSS-REINFORCED WITH STIFFENERS 
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FIGURE 10. STRESS SPECTRAL DENSITY AT CENTER OF TWO OPPOSITE EDGE 
SIMPLY-SUPPORTED AND OTHER TWO CLAMPED RECTANGULAR 
CURVED PANEL CROSS-REINFORCED WITH STIFFENERS 
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FIGURE 11 ACCELERATION SPECTRAL DENSITY AT CENTER OF PANEL, 
ALL TEEMS SUMMATION 
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FIGURE 12 ACCELERATION SPECTRAL DENSITY AT CENTER OF PANEL, 
CROSS-TERMS NEGLECTED 
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FIGURE 13 DISPLACEMENT SPECTRAL DENSITY AT CENTER OF 
PANEL, ALL TERMS SUMMATION 
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FIGURE 14 DISPLACENT SPECTRAL DENSITY AT CENTER OF PANEL, 
CROSS-TERMS NEGLECTED 
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FIGURE 15 STRESS SPECTRAL DENSITY AT CENTER OF PANEL ALL 
TERMS SUMMATION 
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FIGURE 16 STRESS SPECTRAL DENSITY AT CENTER OF PANEL, 
CROSS-TERMS NEGLECTED 
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FIGURE 17 ACCELERATION SPECTRAL DENSITY, ALL TERMS SUMMATION, 

AT x = %! 4 = 11.87 in. y = b/4 = 14.594 in. 
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FIGURE 19 


DISPIACEMENT SPECTRAL DENSITY, ALL TERMS SUMMATION, 
AT x a ft/4 = 11.87 In. y - b/4 = 14.594 in. 
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FIGURE 20 DISPLACEMENT SPECTRAL DENSITY, CROSS-TERMS NEGLECTED 
AT x - A/4 = 11.87 in. y = b/4 = 14.594 in.’ 
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FIGURE 21 STRESS SPECTRAL DENSITY, ALL TERMS SUMMATION, 

AT x = £/4 = 11.87 in. y = b/4 = 14.594 in. 
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FIGURE 22 


STRESS SPECTRAL DENSITY, CROSS-TERMS NEGLECTED, 

AT x - £/4 — 11.87 in. y — b/4 = 14.594 in 
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FIGURE 24 ACCELERATION SPECTRAE DENSITY, CROSS-TERMS NEGLECTED 
AT x = A - 23.75 INCHES, y = A = 14.594 INCHES 
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FIGURE 25 DISPLACEMENT SPECTRAL DENSITY, ALL TERMS SUMMATION, 

AT x = A = 23.75 INCHES, y =.]L = 14.594 INCHES 
2 4 
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FIGURE 26 DISPLACEMENT SPECTRAL DENSITY, CROSS-TERM NEGLECTED, 

AT x = y = 23.75 INCHES, y = * = 14.594 INCHES 

z 4 
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FIGURE 27 STRESS SPECTRAL DENSITY, ALL TERMS SUMMATION, 
AT x = ~ = 23.75 INCHES, y = ~ = 14.594 INGHES 
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FIGURE 29 


ACCELERATION SPECTRAL DENSITY, ALL TERMS SUMMATION, 
AT x = V4 = 11. 87 in . , y = b/2 = 29. 187 in. 
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FIGURE 31 DISPLACEMENT SPECTRAL DENSITY, ALL TERMS SUMMATION, 
AT x - A /4 = 11. 87 in. , y = b/2 = 29. 187 in. 
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FIGURE 33 STRESS SPECTRAL DENSITY, ALL TERMS SUMMATION, 
AT x = V4 = 11. ,87 in. , y = b/2 = 29. 187 in. 
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FIGURE 35 ACCELERATION SPECTRAL DENSITY AT CENTER OF PANEL, 
RADIUS OF CURVATURE a = 100 INCHES 
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FIGURE 37 DISPLACEMENT SPECTRAL DENSITY AT CENTER OF PANEL, 
RADIUS OF CURVATURE a = 100 INCHES 
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FIGURE 38 AVERAGE DISPLACEMENT SPECTRAL DENSITY OVER WHOLE 
PANEL, RADIUS OF CURVATURE a = 100 INCHES 
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TOTAL NUMBER OF MODES 


PANEL DATA: LENGTH = 47. 5 INCHES 
WIDTH = 58. 4 INCHES 
THICKNESS = 0. 75 INCH 
YOUNG'S MODULUS = 10 X 10 6 lbf/in 2 
POISSON'S RATIO = 0.287 



FIGURE 41 NUMBER OF MODES AS A FUNCTION OF FREQUENCY OF A SIMPLY 
SUPPORTED UNIFORM RECTANGULAR PANEL 
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MODAL DENSITY IN MODES PER HERTZ 


PANEL DATA: LENGTH = 47. 5 INCHES 
WIDTH = 58.4 INCHES 
THICKNESS = 0. 75 INCH 
YOUNG’S MODULUS = 10X10 6 lb£/in 2 
POISSON’S RATIO = 0. 287 



FIGURE 42 MODAL DENSITY OF A SIMPLY SUPPORTED 
UNIFORM RECTANGULAR PANEL 
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f = 100 Hertz 



SEPARATION 0 1 STANCE UNChESI 


FIGURE 43. PLOT OF EXPONENTIALLY DECAYING CORRELATION FUNCTION WITH 
DECAY CONSTANT A = 5 
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FIGURE 44. PLOT OF EXPONENTIALLY DECAYING CORRELATION FUNCTION WITH 
DECAY CONSTANT A - 0.5 
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FIGURE 45. PLOT OF SINUSOIDAL DECAYING CORRELATION FUNCTION WITH DECAY 
CONSTANT A = 5 
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c = 13500 in. /sec. 

f = 100 Hertz 



FIGURE 46. PLOT OF SINUSOIDAL DECAYING CORRELATION FUNCTION WITH DECAY 
CONSTANT A' = 0.5 
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FIGURE 47. JOINT ACCEPTANCE, EXPONENTIALLY DECAYING CORRELATION FUNCTION 
MODE INDICES: j=l, k-I, m-1, n=l 
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FIGURE 48. JOINT ACCEPTANCE, EXPONENTIAL DECAYING CORRELATION FUNCTION, 
MODE INDICES: j » 1, k = 3, m = 1, n=3 
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FIGURE 49. JOINT ACCEPTANCE, EXPONENTIALLY DECAYING CORRELATION FUNCTION, 
MODE INDICES: j=l, k=l, m=l, n=3 
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FIGURE 50. NORMALIZED SPECTRAL DENSITY OF GENERALIZED FORGE, EXPONENTIALLY DECAYING 
CORRELATION FUNCTION, MODE INDICES: 1, 1, 1, 3 
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FIGURE 51. JOINT ACCEPTANCE, EXPONENTIALLY DECAYING CORRELATION FUNCTION 
MODE INDICES: i=1k=3.m=1.n=S 
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FIGURE 53. JOINT ACCEPTANCE, SINUSOIDAL DECAYING CORRELATION FUNCTION, 
MODE INDICES: j-I.k-l, m-1, n=l 
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FIGURE 54. JOINT ACCEPTANCE, SINUSOIDAL DECAYING CORRELATION FUNCTION, 
MODE INDICES: j = I, k a 1, m= 1, n “ 3 
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FIGURE 55. NORMALIZED CROSS SPECTRAL DENSITY OF GENERALIZED FORCE, 
SINUSOIDAL DECAYING CORRELATION FUNCTION, MODE INDICES: 
1 , 1 , 1 , 3 
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FIGURE 56. JOINT ACCEPTANCE, SINUSOIDAL DECAYING CORRELATION FUNCTION, 
MODE INDICES: j = 1, k= 3, m = 1, n = 3 
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FIGURE 57. JOINT ACCEPTANCE, SINUSOIDAL DECAYING CORRELATION FUNCTION, 
MODE INDICES: j=l, k - 3, m - 1, n=5 
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FIGURE 58. NORMALIZED GROSS SPECTRAL DENSITY OF GENERALIZED FORCE, 
SINUSOIDAL DECAYING CORRELATION FUNCTION, MODE INDICES: 
1, 3, 1, 5 
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FIGURE 59. DECIBEL SCALE ACCELERATION SPECTRAL DENSITY AT CENTER OF 

SIMPLY-SUPPORTED CURVED RECTANGULAR PANEL CROSS-REINFORCED 
WITH STIFFENERS 
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FIGURE 60' VIBRO-ACOUSTIC TRANSFER FUNCTION AT CENTER OF SIMPLY- 
SUPPORTED CURVED RECTANGULAR PANEL CROSS-REINFORCED 
WITH STIFFENERS 
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FIGURE 61- ACOUSTIC PRESSURE SPECTRAL DENSITY IN DECIBELS FOR FATIGUE 
TEST IN PROJECT NAS8-21425 
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FIGURE 62 . COMPARISON OF CALCULATED RESULT WITH 
TEST RESULT, LOCATION NO. 9 
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